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Abstract. Over the past two decades, unconventional superconductivity with 
gap symmetry other than s-wave has been found in several classes of materials, 
including heavy fermion (HF), high-T c , and organic superconductors. Unconventional 
superconductivity is characterized by anisotropic superconducting gap functions, which 
may have zeros (nodes) along certain directions in the Brillouin zone. The nodal 
structure is closely related to the pairing interaction, and it is widely believed 
that the presence of nodes is a signature of magnetic or some other exotic, rather 
than conventional phonon-mediated, pairing mechanism. Therefore experimental 
determination of the gap function is of fundamental importance. However, the detailed 
gap structure, especially the direction of the nodes, is an unresolved issue in most 
unconventional superconductors. Recently it has been demonstrated that the thermal 
conductivity and specific heat measurements under magnetic field rotated relative 
to the crystal axes are a powerful method for determining the shape of the gap 
and the nodal directions in the bulk. Here we review the theoretical underpinnings 
of the method and the results for the nodal structure of several unconventional 
superconductors, including borocarbide YM2B2C, heavy fermions UPd2Al3, CeCoIns, 
and PrOs4Sbi2, organic superconductor, k-(BEDT-TTF) 2 Cu(NCS)2, and ruthenate 
Si'2Ru04, determined by angular variation of the thermal conductivity and heat 
capacity. 
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1. Introduction 

Superconducting transition is a second order phase transition associated with a 
spontaneous symmetry breaking. Consequently the order parameter that appears below 
the transition temperature, T c , characterizes the lowering of the symmetry in the low 
temperature ordered phase compared to the metallic state. The order parameter is also 
related to the gap in the single particle excitation spectrum, and hence its symmetry 
is reflected in the elementary excitations in the superconducting phase. These, in turn, 
determine the observed transport and thermodynamic properties jTJ [2] . 

At the microscopic level, the symmetry of the order parameter is intimately related 
to the pairing interaction. Therefore the identification of this symmetry is of central 
importance in the study of the superconductivity, especially in novel correlated electron 
materials. This review describes recent progress and current status of the efforts to use 
transport properties as a reliable tool for the determination of the gap symmetry in 
unconventional superconductors in the bulk. 

In all superconductors the gauge symmetry is broken below the transition 
temperature. At the same time, in most materials, the energy gap has the full spatial 
symmetry of the underlying crystal lattice. In the simplest cases, this corresponds to a 
gap isotropic in the momentum space, i.e. independent of the directions at the Fermi 
surface. These superconductors are termed conventional, or s-wave. However, it is also 
possible that the spatial symmetry of the superconducting order parameter is lower 
than that of the lattice, and such superconductors are labeled unconventional. Secj2] 
presents the detailed symmetry classification of the unconventional superconductors. 
These materials first appeared in 1979, when the heavy fermion CeCu2Si2 was discovered 
[3], and by now superconductivity with non s-wave symmetry is ubiquitous. Examples of 
it include anisotropic gaps with zeroes (nodes), odd parity superconducting condensate 
wave functions, and broken time reversal symmetry. Realization of these possibilities 
is proved or strongly suggested in several classes of materials, including heavy 
fermion [USEE], high-T c cuprates [7], ruthenate[8| 191 HO], cobaltate [TT], intermetallic 
compounds [12], and organic superconductors [T3| [HI [151 [16] . 

In these systems, strong electron correlations often give rise to the Cooper 
pair states with a non-zero angular momentum. Unfortunately, the experimental 
determination of the detailed superconducting gap structure is an extremely difficult 
task. The phase sensitive measurements testing the sign change of the gap have only 
been done in the high-T c cuprates, firmly establishing (along with a number of other 
techniques probing the gap anisotropy) the predominant d x i_ y i pairing symmetry [7]. 

Unconventional superconductivity often occurs in heavy fermion compounds, 
containing /-electrons (lanthanide (4f) and actinide (5f)), especially in materials 
containing Ce, Pr, U and Pu atoms. At high temperature /-electrons are essentially 
localized with well-defined magnetic moments. As the temperature is lowered, the 
/-electrons begin to delocalize due to the hybridization with conduction electron 
band(s,p, (i-orbital), and Kondo screening. At yet lower T, the /-electrons become 
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itinerant, forming a narrow conduction band with heavy effective electron mass (up to 
a few hundred to a thousand times the free electron mass). Strong Coulomb repulsion 
within a narrow band and the magnetic interaction between remnant unscreened 4/ 
or 5/ moments leads to notable many-body effects, and, likely, to superconductivity 
mediated by magnetic fluctuations. These effects are believed to be especially 
pronounced in the vicinity of zero-temperature magnetic instability (quantum critical 
point (QCP)). Although the superconducting gap is thought to be unconventional 
in most heavy fermion superconductors, its detailed momentum dependence is an 
unresolved issue [U [5j [6] . 

Since the discovery of superconductivity in organic materials about two decades 
ago, superconductivity has been reported in more than 100 organic compounds. Among 
them, two families of compounds, quasi-lD Bechgaard salts (TMTSF) 2 X (X=C104, PFq, 
AsFg, etc.) and quasi-2D k-(BEDT-TTF) 2 X salts (/c-(ET)2X), where the ion X can, for 
example, be Cu(SCN) 2 , Cu[N(CN) 2 ]Br or I 3 , are particularly prominent candidates for 
unconventional pairing. The pairing symmetry in both systems is still undetermined, 
and is one of the most intriguing problems in the field [HI ISl EE!] • 

The unconventional superconductivity has also been reported in some of transition 
metal oxides other than high-T c cuprates. Especially, the superconducting gap functions 
in layered ruthenate Sr 2 Ru0 4 [9j [10] and layered cobaltate Na x Co0 2 -yH 2 0[lT] have 
attracted considerable interest. Moreover, among intermetallic compound, the gap 
function of borocarbide superconductors (Y, Lu)Ni 2 B 2 C [12] has been reported to be 
very anisotropic, which implies that the simple electron-phonon pairing mechanism 
originally envisaged for these compounds is not the sole source for pairing interaction. 

Evidence for anisotropic gap in a variety of systems has continued to motivate 
theorists to propose new models for the unconventional superconductivity, which make 
specific predictions for the shape of the superconducting gap in momentum space. 
Experimental determination of the gap symmetry is therefore of crucial importance. 
The transport measurements are not, per se, phase sensitive, and therefore cannot 
unequivocally determine the sign change in the gap function. They are, however, an 
extremely sensitive probe of the anisotropy of the gap amplitude in the momentum 
space, and have been extensively used in the last few years to determine the shape 
of the gap in many materials at the forefront of modern research. Below we give an 
overview of these efforts. 

2. unconventional pairing state 

The general classification scheme for the superconducting order parameter is based on 
its behavior under symmetry transformations. The full symmetry group Q of the crystal 
contains the gauge group U(l), crystal point group G, spin rotation group £77(2), and 
time reversal symmetry group T, 

Q = U{l)®G®SU{2)®T . (1) 
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The superconducting order breaks the U(l) gauge symmetry below T c , and the simplest 
superconductors are those in which only the U(l) symmetry is broken; these are labeled 
conventional. Unconventional superconductors break an additional symmetry besides 
U(l), and may include order parameters which 

(i) have odd parity; 

(ii) break time reversal symmetry; 

(iii) break the point group symmetry of the crystal. 

The superconducting order parameter is proportional to the gap function Af S2 (k), 
which, in turn, is proportional to the amplitude of the wave function for a Cooper 
pair \Pg (k) = {ipk^ip -k,s 2 )- Here k is the quasiparticle momentum, Sj is the 
electron spin, and is the electron annihilation operator. The order parameter is 
called unconventional if it transforms according to a nontrivial representation of the 
full symmetry group. Pauli's exclusion principle requires Ag S2 (k) to be antisymmetric 
under the the particle interchange: Af S2 (k) = — Ag 2 Si (— k). In the simplest case of 
weak spin-orbit interaction, the total angular momentum L and total spin S = Si + s 2 
are good quantum numbers, and A Sl;S2 (k) can be written as a product of orbital and 
spin parts, 



The orbital part, g^(k), can be expanded in the spherical harmonics Y^ m (k), which are 
the eigenfunctions of the angular momentum operator with the momentum £ and its 
^-projections m, 



Here k = k/fcp denotes the direction of the Fermi surface. ge(k) is even for even 
values of £ and odd for odd values of £, ge(k) = (— k), and superconductors 
with £ = 0,1,2,3,4,... are labeled as having s,p,d, f, g, . . . wave gap respectively. 
This classification is valid for an isotropic system. In a crystal, the spatial part of 
the Cooper pair wave function is classified according to the irreducible representations 
of the symmetry group of the lattice. However it is common even in this case to refer 
to the possible pairing states as having a particular angular momentum (rather than 
belonging to a representation of the group with given symmetry properties) and we use 
the notation here. 

The spin part of the order parameter, Xs( s i, s 2), is a product of the spinors for the 
two electrons in the Cooper pair. Therefore the gap function is a 2x2 matrix in spin 

space, 




(2) 




(3) 



m=-£ 





5 



In the singlet state, S = 0, the spin part of the wave function is | t!) — I !t)> an d 
therefore the gap function is simply proportional to the Pauli matrix a y : 

A s (k) = Ag e (k)ia y , (5) 

where £ is even and gg is normalized. The energy of single particle excitations in this 



case is 



£ k = ^ + A2 |^ (k) | 2; (6) 

where £k is the band energy relative to the chemical potential. For superconductors 
with an isotropic A(k) the excitations have a finite energy gap everywhere at the Fermi 
surface, while for anisotropic pairing the gap amplitude depends on the components of 
<?(k). 

For spin triplet pairing (S = 1), the wave function has components corresponding to 
the three different spin projections, S z , on the quantization axis (|Tt)) It!) + |!t); I!!))- 
Consequently, it is common to write the order parameter as 

A t (k) = iA(d(k)a)a y 

( -d x (k) + id y (k) 4(k) \ ( 

\ 4(k) d X (k)+ldy(k) J ^ ) 

The orbital part is expressed by these d-vector with 

gi = ~d x + id y , g 2 = d z , g 3 = d x + id y , (8) 
and the excitation energy is 



£k = ^ + A 2 |d(k)p. (9) 

In the presence of strong spin-orbit coupling only the total angular momentum 
J=L+S is a good quantum number, and the classification according to physical electron 
spin is not possible. However, if the crystal structure has inversion center, Cooper 
pair states can still be classified according to their parity, and therefore acquire a 
"pseudospin" quantum number, instead of the physical spin. This situation is commonly 
encountered in many heavy fermion materials. 

Experimentally, the parity of the pair wave function ty pair (k) can be determined 
by the Knight shift of the nuclear magnetic resonance (NMR) frequency [2j[T7], muon 
spin rotation (/uSR), and, less directly, by the magnitude of the upper critical field 
H C 2. The Knight shift is linear in the electron spin susceptibility Xs, and is therefore a 
direct measure of the spin polarization in the superconducting state. In a spin singlet 
superconductor, as Cooper pairs are formed, the spin contribution to the Knight shift 
falls rapidly on cooling through the transition. In contrast, in a triplet superconductor 
the spin orientation of the Cooper pairs is determined by the d-vector in Eq. ([TJ) . If the 
direction d is fixed by the spin-orbit interaction, the Knight shift is anisotropic in the 
superconducting state. When the magnetic field is applied along d (H|| d), the Cooper 
pair spin is perpendicular to H and hence does not contribute to the susceptibility. 
Then the Knight shift decreases rapidly below T c , as in spin singlet superconductors. 
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On the other hand, when the applied field H_L d, the x- and y- components of d vector, 
d x (k) + id y (k), are nonzero, and the contribution to the susceptibility of the Cooper 
pairs is identical to that of the constituent electrons, x\\ — Xn, where n labels the normal 
state above the transition temperature. Therefore the Knight shift remains unchanged 
below T c . 

Applied magnetic field destroys superconductivity through both the orbital 
dephasing and Zeeman splitting of the single electron energy levels. In type-II 
superconductors the former effect leads to the emergence of the vortex state, when the 
magnetic field penetrates into the sample, forming the regular array of the vortex tubes 
parallel to the field. Vortices have cores of the size of the superconducting coherence 
length, £, where superconductivity is destroyed. Each vortex carries a flux quantum 
$o — Tthc/e = 2 x 10 G-cm 2 , and therefore in the external field, H, the area per vortex, 
A, is determined from AH = $o- At the orbital upper critical field the vortex cores 
overlap destroying bulk superconductivity, hence a simple estimate gives H° 2 = $o/2vr^ 2 . 
On the other hand, in singlet superconductors, an additional pairbreaking effect of 
the field is due to polarization of the normal state electrons. The upper critical field 
determined by this Pauli limiting effect is estimated to be H^ 2 = A/\/2/xb, where /xg is 
the Bohr magneton, and A is superconducting gap. The Pauli limiting is absent in spin 
triplet superconductors. Therefore, finding H c2 which is higher than H^ 2 may indicate 
spin triplet pairing. 

Up to now, possible odd parity superconducting state has been suggested in 
heavy fermion UPt 3 [18j, UNi 2 Al 3 [19j, URu 2 Si 2 |2Q|, UBei 3 ^T], PrOs 4 Sbi 2 [25], organic 
(TMTSF) 2 PF 6 |23j, and transition metal oxides, Sr 2 Ru0 4 [2I], Sr 2 Cai 2 Cu 24 4 i[25] and 
Na x Co0 2 -yH 2 [26j[27]. It is most probably realized in ferromagnetic superconductors 
UGe 2 [28], URhGe[29] and UIrj30]. However, one needs to bear in mind that the odd 
parity in some of these materials is still controversial. In fact, analysis of the NMR 
spectrum in the vortex state, from which we determine the Knight shift, is not settled. 
Moreover the NMR experiment measures the surface area with the length scale of 
penetration depth A . In this regime, strong currents associated with the surface barrier 
flow and the field distribution and the susceptibility is strongly space inhomogeneous. 
The influence of these currents on the NMR spectrum is an open question. In addition, 
if the pairing interaction is modified by magnetic field, the upper critical field can be 
enhanced above H^ 2 even in the spin singlet superconductors. 

When typ a i r (k) has an imaginary part 

^ r (k) =Vi(k)+^ 2 (k), (10) 

the time reversal symmetry is broken, since ^ Pair (k) ^ ^ r pa j r (k). In such a situation, 
spontaneous static magnetic field arising from the orbital current around the impurity 
or at the surface can appear below T c , because the impurity or boundary lifts the 
degeneracy between \L,L Z ) and \L, — L z ). Such a spontaneous magnetic field was 
observed in UPt 3 [31], Sr 2 Ru0 4 |32j, and PrOs 4 Sb 12 [33] by /iSR experiments. 
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Recently superconductivity with no spatial inversion symmetry has excited great 
interest. In the presence of strong spin-orbit interaction, the absence of the inversion 
symmetry strongly influence the pairing symmetry through a splitting of the two spin 
degenerate bands. Generally in the system without inversion symmetry the gap function 
is a mixture of spin singlet and triplet channels in the presence of a finite spin orbit 
coupling strength. Associated with the absence of spatial inversion symmetry, unusual 
superconduting properties, including the striking enhancement of H c2 and helical vortex 
phase, have been proposed. The absence of inversion symmetry has been reported in 
several superconductors. Among them, the superconducting gap function with line 
nodes have been reported in in CePt3Si[34j and I^PtsB [35]. In these system, the 
position of node is strongly influenced by the mixing of the spin singlet and triplet 
components. 

In Table-I, we summarize the gap functions of several anisotropic superconductors. 
One needs to bear in mind that the gap functions in some of the listed materials are 
still controversial. 

3. Nodal Structure: Standard Techniques 

In unconventional superconductors discovered so far, the energy gap for the quasiparticle 
excitations vanishes (has nodes) at points or along lines on the Fermi surface. There is 
now a wide variety of dynamic and thermodynamic probes that couple to, and reveal 
these low energy excitations. The temperature dependence of the London penetration 
depth A(T), electronic part of the specific heat C(T), thermal conductivity n(T), 
and nuclear magnetic resonance (NMR) spin-lattice relaxation rate T£~ all reflect 
the changes in the quasiparticle occupation numbers. In the fully gapped (s-wave) 
superconductors the quasiparticle density of states (DOS) is zero at energies, E, below 
the gap edge (no excitations with E < A), and varies as N S (E)/N = y / E i_ A i for E > A. 
Here Nq is the normal state DOS. The physical quantities exhibit activated temperature 
dependence, exp(— A/T) at low temperatures, T <C T c . On the other hand, in nodal 
superconductors, the low-energy density of states remains finite due to contributions 
from the near-nodal regions on the Fermi surface, and typically the DOS varies as 
N S (E)/N = (E/A ) n at E C A. The exponent n depends on the topology of the 
nodes: n = 1 for line nodes as well as for point nodes where the gap is quadratic 
in distance from the nodal point in the momentum space; n = 2 for point nodes 
with linearly varying gap amplitude around the nodal point. Then the experimental 
quantities described above exhibit power law temperature dependence at T T c . For 
example, in <i-wave superconductors with line nodes, the DOS N(E) ~ \E\ leads to 
the specific heat C e ~ 7 n T 2 /T c , where 7 n is the coefficient of the linear-T-term in the 
normal state, and the NMR relaxation rate T^ 1 oc T 3 . The deviation of the superfluid 
density, n s (T) from its zero temperature value, An s (T) = n s (0) — n s {T) oc ksT/Ao, 
which can be detected by the penetration depth measurements. 

So far we discussed pure systems. The regime where power laws in T are 
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observed is strongly influenced by the impurities. In unconventional superconductors 
non-magnetic impurities act as pair-breakers, similar to magnetic impurities in s- 
wave superconductors. A bound state appears near an isolated non-magnetic strong 
(scattering phase shift tt/2, or unitarity) scatterer, at the energy close to the Fermi 
level. The broadening of this bound state to an impurity band at finite disorder leads 
to a finite density of states at zero energy, N(0), that increases with increasing impurity 
concentration [36] . The impurity scattering changes the temperature dependence of the 
physical quantities below T corresponding to the impurity bandwidth: AA changes 
the behavior from T to T 2 , the NMR relaxation rate changes from T x ~ oc T 3 to 
T-j -1 oc T (T]T ~ const.), and C(T) changes from T 2 to T. In NMR, for example, 
the temperature range where T^ 1 exhibits the T 3 - dependence is limited to T > T c /3 
in most measurements. 

Therefore the nodal behavior may be hidden by the impurity effects. Even in the 
extremely pure systems, however, experimental observation of the power laws provides 
indications of the existence of the nodes, but is unable to yield information about their 
location in the momentum space. 

Angle resolved photoemission spectroscopy (ARPES) directly investigates the 
momentum dependence of the gap, and was instrumental in determining the gap shape in 
high-T c superconductors. However, the energy resolution of this technique is insufficient 
when compared with the size of the energy gap in most low-T c systems. The phase- 
sensitive measurements, such as corner junctions, tricrystal, and tests for Andreev bound 
states, which provided the most convincing evidence for d x 2_ y 2 wave order parameter 
in the high-T c cuprates, are primarily surface probes. Absence of inversion symmetry 
near the surface may influence the pairing symmetry through, for example, splitting 
of the two spin degenerate bands via the spin-orbit coupling, so that the gap function 
is a mixture of the spin-singlet and spin-triplet channels. Moreover, it is difficult to 
apply these techniques to determine the three dimensional gap structure, and they have 
received limited use beyond studies of the high-T c cuprates. It is therefore extremely 
important to acquire complementary evidence for particular gap symmetries via bulk 
measurements. 

4. Nodal Superconductor in Magnetic Field 

4-1. General approaches to the vortex state. 

Determining the nodal positions requires a directional probe. In the following we argue 
that an applied magnetic field provides a convenient bulk probe of the symmetry of the 
order parameter in unconventional superconductors. The usefulness of the magnetic 
field as a probe relies on an important difference between the properties of the vortex 
state in nodal compared to fully gapped s-wave superconductors. While for the s-wave 
case the DOS and the entropy at low fields, H <C H c2 , are determined by the localized 
states in the vortex cores, in the superconductors with nodes they are dominated by 
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the extended quasiparticle states, which exist in the bulk close to the nodal directions 
in momentum space. Therefore much attention has been paid to the effect of the field 
on these near-nodal quasiparticle. 

A simple picture that captures the main effect of the magnetic field on a nodal 
superconductor is that of the Doppler shift of the quasiparticle spectrum [37J. In the 
presence of a supercurrent with velocity v s the energy of a quasiparticle with momentum 
k is Doppler shifted relative to the superconducting condensate by 



This effect originates from the Galilean transformation: creation of an excitation (p, e) 
in the normal-quasiparticle rest frame, involves an additional energy 5e = p ■ v s in the 
superfluid frame of reference. For a uniform superflow the Doppler shift is simply a 
consequence of the gauge invariance, and is therefore exact [38]. For a non-uniform 
superflow, as in the vortex state, this picture is semiclassical in that it considers 
simultaneously the momentum of the quasiparticles and the local value of v s (r) at 
position r, and therefore ignores the possible accumulation of the quantum mechanical 
phase around the magnetic vortices in superconductors. The fully quantum mechanical 
treatment of the quasiparticle energies so far was carried out only in a perfectly periodic 
vortex lattice and in the absence of impurities [39], and gives results for the physical 
properties close to those obtained in the semiclassical treatment [391 SO] • 

To estimate the characteristic energy scale of the Doppler shift we can approximate 
the velocity field by that around a single vortex, v s = h(f)/2mr, where r is the distance 
from the center of the vortex and is a unit vector along the circulating current. This 
expression is valid outside the vortex core and up to a cutoff of order min{i?, A}, where 
R = ay/Qo/nH is the intervortex distance, <3>o 1S the flux quantm, a is a geometric 
constant, and A is the London penetration depth. Average Doppler shift, E av , is 
computed by integrating over a vortex lattice unit cell, and is given by 



Thus E av is proportional to 

Since the density of states is an additive quantity, the net DOS of the sample is 
the sum of the contributions from the areas with distinct values of the Doppler shift. 
In a system with line nodes, where the low energy DOS N(E) oc \E\, this implies the 
residual density of states N V (E = 0,H) oc E av oc yH. Consequently the specific heat 
also exhibits the V^ff-behavior at low temperatures in the clean limit [37J HU H21 H3]. 
Since the supervelocity distribution can be obtained for a given configuration of vortices, 
the range of possible values for the coefficient of the \/~H can also be found for a given 
material [44J. 

Determining the transport properties, such as the thermal conductivity k, using 
the Doppler shift method is a more challenging task. The transport coefficients are 
determined from the correlation functions that have a finite range, and therefore depend 
on the Doppler shift at more than a single point. Local values of these coefficients can be 



e(k) -> e(k) -frk-v. 



(11) 




(12) 
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rigorously defined only in the dirty limit. It is generally accepted that a similar definition 
gives at least a qualitatively correct results in the clean limit [15], although a rigorous 
comparison is currently lacking. Even with that assumption, the connection between 
the distribution of local values, k(t), and the measured value remains a subject of some 
debate. Both averaging k(t) and « _1 (r) have been proposed [H2, 06]. While in some 
cases the difference is only in the magnitude of the field-induced change, the divergent 
philosophies behind the averaging procedures give rise to qualitatively different results 
for the anisotropy of the transport coefficients described below. 

In the approach discussed above only the quasiparticle energy is shifted, so that 
the single particle scattering rate is not directly affected by the presence of the vortices. 
In the presence of static disorder treated, for example, in the self-consistent T-matrix 
approximation, the magnetic field does affect the lifetime indirectly, by modifying the 
density of states available for scattering [HJ 05] . Hence the Doppler shift method does 
not account for scattering of the quasiparticles on vortices. 

Calculations of the vortex scattering cross-section have to go beyond the 
semiclassical treatment and make assumptions about the structure of the vortex core 
states [47] . and therefore received limited attention. An alternative, fully microscopic 
approach, employs an extension of the approximation originally due to Brandt, Pesch 
and Tewordt (BPT) [HJ HH] to describe clean superconductors near the upper critical 
field, H C 2. In this method the Gor'kov equations (or their quasiclassical Eilenberger- 
Larkin-Ovchinnikov analog) are solved with the normal electron Green's function 
replaced by its spatial average [33], EH EU E21 [53]. This method is rigorously justified 
at moderate to high fields, gives the standard quasiparticle spectrum as H — > 0, and 
yields results that are qualitatively similar to the Doppler shift at very low fields and 
temperatures. Therefore it is believed that it can be used over a wide range of fields 
relevant to experiment. The BPT approach naturally includes the scattering of the 
quasiparticles off the vortices [50, [53]; however, due to the incoherent averaging over 
different unit cells of the vortex lattice, it tends to overestimate the importance of 
such scattering at lower fields. Together, BPT and Doppler shift methods account 
for a majority of theoretical work relevant to the experimental investigations of the 
quasiparticle properties in the vortex state of nodal superconductors. 

4-2. Thermal Conductivity 

Of all the transport properties the thermal conductivity is uniquely suitable for 
probing bulk superconductivity. Unlike electrical resistivity, it does not vanish in 
the superconducting state. Cooper pairs do not carry entropy and therefore do not 
contribute to the thermal transport. As a result, the thermal conductivity probes the 
delocalized low energy quasiparticle excitations, and is sensitive to the effect of magnetic 
field on the quasiparticles. In Fig. [1] we show the qualitative behavior of the thermal 
conductivity and heat capacity at low T as a function of the magnetic field for a s-wave 
(fully gapped) and a d-wave (with line nodes) superconductor. 



12 




Figure 1. Field dependence of (a) the thermal conductivity k and (b) the specific 
heat C for s- and d-wave superconductors. In s-wave superconductors, the thermal 
conductivity shows an exponential behavior with very slow growth with H [541 155] . 
The heat capacity increases nearly linearly with H. In sharp contrast, in ci-wave 
superconductors, both the specific heat and the quasiparticle conduction grows rapidly 
as soon as the field exceeds H c %. The slope of k(H) at H C 2 depends on purity and 
therefore for the (i-wave case it is also possible to have an inflection point in k(H) at 
intermediate fields. 



In s-wave superconductors the only quasiparticle states present at T C r c are 
those associated with vortices. At low fields where the vortices are far apart, these 
states are bound to the vortex core and are therefore localized and unable to transport 
heat; the thermal conductivity shows an exponential behavior with very slow growth 
with H. At high fields near H C 2 where quasiparticles states within the vortices begins 
to overlap with those within the neighboring vortices, thermal conductivity increases 
rapidly. Such a field dependence of the thermal conductivity is observed in Nb [54l 155]. 
The heat capacity, due to the localized quasiparticle states, increases nearly linearly 
with H. In dramatic contrast, both the specific heat and the quasiparticle conduction, 
due to near nodal states, grow rapidly as soon as the field exceeds H c \. In <i-wave 
superconductors where N(E = 0, H) oc y/H due to the Doppler shift of the quasiparticle 
energy spectrum, both the thermal conductivity and the specific heat exhibit a nearly 
v^ff-behavior. 

In reality, especially at higher temperatures, the behavior of the thermal 
conductivity is more complex. While the magnetic field enhances the local DOS, 
it also leads to a change in the transport lifetime both via the modification of the 
impurity scattering and via Andreev scattering off the vortices. Understanding of these 
competing effects has progressed during the past few years [HI HH [50l [53j EH [57J [581 [59] , 
although the complete picture is not yet developed. In general, at low temperatures the 
DOS modification plays the dominant role, and the thermal conductivity increases with 
increased field. At higher temperatures and low fields, the dominant effect of vortices 
is to introduce and additional scattering mechanism, while the DOS is controlled by T. 
Consequently, the thermal conductivity initially decreases with field, and goes through 
a minimum at a finite H [50]. This behavior has been first observed in high-T c cuprates 
[60] , and also seen in other systems [611 [62] 
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5. How to determine the nodal structure in the bulk? 

5.1. Anisotropy under rotating field: density of states 

The techniques described above help determine the general topology of the gap, but 
cannot establish the exact angular dependence of A(k). In particular, the nodal 
positions in k-space cannot be obtained. In the following we discuss the theoretical 
underpinnings and experimental realizations of the new and powerful method for 
determining the nodal directions. 

The method is based on the prediction that, under an applied magnetic field, the 
density of states of nodal superconductors depends on the orientation of the field with 
respect to the nodal direction [63], and exhibits characteristic oscillations as the field 
is rotated relative to the sample. The oscillations can be measured via the field-angle 
dependence of the thermal conductivity [58[ EH EH [62], E3 E3 E3 EH EH] or specific 

heat isa cna eh E2i csi . 

While the specific heat measurements directly probe the density of states, the 
thermal conductivity anisotropy is sensitive to a combination of the density of states, 
and the quasiparticle transport scattering rate, which may have different dependence on 
the field orientation. Measurements of the specific heat anisotropy were only attempted 
several years after theoretical predictions [69], [70], EH E21 [73]. First experiments on 
the field-angle dependence of the thermal conductivity preceded theoretical discussions 
[581 [59], but focused simply on the existence, rather than location, of additional features 
(interpreted as arising from Andreev scattering) in cuprates. Use of the thermal 
conductivity as a similar test based on the DOS anisotropy was, to our knowledge, first 
suggested in Ref. [71], and followed up by other work [16]. Development and consistent 
use of the measurements to probe the direction and type of nodes in k-space, is almost 
entirely due to recent efforts by the group of University of Tokyo [HH E2[ [65], EHl [67] [68] . 
The full theory of the anisotropy of the thermal conductivity in the vortex state is still 
incomplete. Hence, while the salient features of experiments are qualitatively understood 
based on a number of treatments [16] [53] [63] [75] [76] [77] [78] [79] [80] EH E21 ES], many 
details need to be addressed further. Below we discuss the current status of this field. 

The origin of the anisotropy is best understood in the framework of the Doppler shift 
of the delocalized quasiparticle spectrum in the vortex state. Consider, for simplicity, 
a d xy gap symmetry with four vertical lines of nodes, and assume a cylindrical or 
spherical Fermi surface, as illustrated in FigsfS] (a)-(c). At low fields, the loci of 
unpaired quasiparticles in the momentum space are close to the nodal lines. Since 
the supercurrents flow in the plane normal to direction of the applied field, the Doppler 
shift experienced by quasiparticles in a given near-nodal region depends on the direction 
of the field H(9, 0) = if (sin 9 cos 0, sin 9 sin 0, cos 9) with respect to the nodal directions. 
Here 9 and are the polar angle and the azimuthal angle respectively, measured relative 
to the c-axis. 

Consider H rotated conically (fixed 9) with varying in-plane angle, 0; see the view 
from above in Fig. [2(b). When the field is aligned with a nodal line, the superflow 
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Figure 2. (a)Sketch of the gap structure with four line nodes perpendicular to the 
basal plane (vertical node), (b) Schematic diagram showing the regions on the Fermi 
surface that experience the Doppler shift in H within the basal plane. We have 
assumed d xy symmetry. <p = (H,a) is the azimuthal angle measured from the a- 
axis. With H applied along the antinodal directions, all four nodes contribute to the 
DOS, while for H applied parallel to the node directions, the Doppler shift vanishes 
at two of the nodes, (c) Four-fold oscillation of the DOS for H rotating in the basal 
plane. The DOS shows a maximum (minimum) when H is applied in the anitinodal 
(nodal) direction. 

around the vortices is in the plane nearly normal to the momenta of quasiparticles close 
to that node. As a result, for these quasiparticles the Doppler shift is small. In contrast, 
when the field is in the antinodal direction, the Doppler shift is (relatively) large along 
all four nodal lines. As a result, the net DOS has minima when H is aligned with the 
nodal direction, and maxima for H along the antinodes [63J. The angle- dependence of 
the DOS exhibits characteristic four- fold oscillations, as shown in Fig. EJ^c). In general, 
DOS oscillates with n-fold symmetry corresponding to the number of vertical nodes n. 

In this approach the amplitude of the DOS oscillations, 5N(E)/N(E), depends on 
the shape of the Fermi surface and other parameters of the models. For the residual 
(E = 0) DOS, most calculations predict the oscillation amplitude ranging from 3% to 
10%. The anisotropy is rapidly washed away at finite energy, and therefore the amplitude 
of the corresponding oscillations in the measured quantities, such as the specific heat, 
at finite temperature, is typically of the order of a few percent. 

Consider now horizontal line nodes in a cylindrical or spherical Fermi surface, as 
illustrated in Figj3] (a). The density of states is anisotropic under the rotation of the 
field in the ac-plane, by varying the angle 9. To illustrate the difference between the 
line nodes at high symmetry positions in the Brillouin zone, and away from those, we 
consider here two model gap functions 

(i) type-I : Horizontal nodes located at the center of the Brillouin zone and at the zone 
boundary A(k) oc sin k z c. 

(ii) type-II : Horizontal node located at positions shifted off the zone center; A(k) oc 
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Figure 3. (a) Schematic figure of the gap structure with line nodes parallel to the 
basal plane (horizontal node) in spherical and open Fermi surfaces. Type-I assumes a 
gap function A(k) oc sinA^c, where line nodes are located at the center of the Brillouin 
zone and at the zone boundary. Type-II assumes a gap function A(k) oc cos k z c, where 
line nodes are located at positions shifted off the zone center, (b) Oscillations of the 
DOS for H rotating in the ac-plane for various gap functions. Two-fold oscillation with 
the same sign are expected for type-I. On the other hand, for type-II, an oscillation 
with a double minimum is expected. 

cos k z c. 

The expected angular variation of the Doppler shifted DOS is a function of the relative 
angle between H and p for these gap functions is shown schematically in Figj3] (b). The 
twofold oscillation is expected for type-I gap functions, in which the horizontal nodes 
are located at the position where p|| afr-plane. On the other hand, for type-II, one 
expects an oscillation with a double minimum structure as a function of 6. Note that 
we sketched the DOS for a fixed H/H c2 ; if a measurement is done at a fixed H, the 
anisotropy of H c2 in a quasi-2D systems superimposes an additional two fold component 
on the oscillations, so that for type-II gap the central maximum is distinct from the other 
two. 

The Doppler shift method does not account for the scattering of electrons on 
vortices. Fully microscopic analyses indicate that inclusion of such scattering further 
reduces the amplitude of the DOS anisotropy [Mj- Furthermore, it has been shown very 
recently that at moderate to high fields and temperatures the vortex scattering leads 
to the inversion of the anisotropy: the density of states is greater for the field along 
the nodal directions than for the field along the gap maxima [53] . While this makes 
the analysis of the specific heat data more complicated, it affects the conclusions drawn 
from the analysis of the transport properties less dramatically [53] , as we discuss in next 
section. 

5.2. Anisotropy under rotating field: thermal conductivity 

Our focus in this review is on the determination of the nodal structure via the thermal 
conductivity measurements. The anisotropy of transport coefficients is given by a 
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combined effect of the angular variations of the density of states, and the angle- 
dependent scattering. The latter effect is not yet fully understood. Since in self- 
consistent treatments the scattering rate of quasiparticles off impurities depends on 
the density of states, the impurity scattering under a rotating field acquires the same re- 
fold anisotropy as the DOS. However, depending on the strength of impurity scattering, 
temperature, and the field, the lifetime may exhibit either maxima or minima for the 
field aligned with the nodes. Consequently, the possibility of the inversion of the maxima 
and minima in the T — H superconducting phase diagram (which was unexpected in 
the behavior of the DOS and the specific heat) was anticipated in the anisotropy of the 
transport coefficients. 

It is believed that the field dependence of the thermal conductivity indicates 
whether the lifetime or the density of states effects dominate. In the regime where 
k(H) decreases with increasing field due to field-enhanced scattering, the maxima of the 
anisotropic conductivity are likely to correspond to nodal direction. In contrast, when 
k(H) increases with field, the density of states effects dominate, and the minima of the 
n-fold pattern indicate the nodes. This conjecture is not rigorous [53], but qualitatively 
correct and provides guidance in the situations when no results of microscopic theory are 
available for a given compound. Moreover, in experiment the angle-induced anisotropy 
of k{H) changes sign close to the point where the field dependence has a minimum, 
supporting this view. 

When heat current, jh, is applied in the basal plane, the angle between j h and H 
is varied as the field is rotated. Consequently, the dominant anisotropy observed in 
experiment is that between the transport along and normal to the vortices, i.e. twofold 
[84] . The nodal contribution appears as a smaller effect on this background, as was first 
seen in the high-T c cuprate YBa 2 Cu307_,5 [5H [59], El]- Note that with few exceptions 
[85] Doppler shift does not describe the combined twofold and nodal anisotropy. 

More sophisticated approaches based on the BPT theory give correct shapes of 
the k(4>) curves, and account for most of the observed features. The details of the 
competition between the twofold and the fourfold oscillations depend on the shape 
of the Fermi surface, role of Zeeman splitting, impurity strength and concentration 
etc. Therefore any semi- quantitative comparison of theory and experiment requires 
knowledge of these as an input, and has only been done for few systems. At the same 
time qualitative conclusions about the shape of the gap can still be drawn from the 
simplified analysis, and we review those for the specific compounds discussed below. 

For relatively three-dimensional systems, the current can be applied along the c- 
axis, and the field rotated conically, varying the azimuthal angle 0, and keeping the 
polar angle 9 constant. In that case the relative orientation of the heat current and 
the field remains unchanged, and the oscillations reflect solely the nodal structure. 
Vortex scattering still modifies the amplitude and the sign of these oscillations, but 
the interpretation is greatly simplified by the absence of the dominant twofold term. 

In this geometry it has been predicted that the ^-dependence of the shape and the 
amplitude of the periodic oscillations provide direct information on the type of nodes, 
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Figure 4. The thermal conductivity k zz (4>,9) (the heat current q || z) when the 
magnetic field is rotated conically as a function of <fr, keeping 9 constant, for two 
different types of nodes, (a)point and (b) line. K zz (<p,9) is normalized by k zz (4:5 ,9). 
The corresponding gap functions are illustrated in the insets, (c) Definition of 9 (polar 
angle) and <f> (azimuthal angle). After Ref.[87]- Note that for vertical line nodes the 
relative amplitude of oscillations for 9 between 45° and 90° depends on the shape of 
the Fermi surface, not accounted for in Ref.thalmaki 

point or line [67]. In Figs. H](a) and (b), we compare the angular variation of the thermal 
conductivity k zz (the heat current q || z) when the magnetic field is rotated conically 
as a function of 0, keeping 9 constant, for two different types of nodes calculated from 
the Doppler-shifted QP spectrum, in accordance with Ref.j67J. Here we adopted gap 
functions A(k) = A o sin(20) (cf-wave) for line node, and A(k) = |A {1 — sin 4 9 cos(4</>)} 
for point node. The latter was proposed in Ref . [86] , but is probably not realized in this 
system, and we use it as a convenient ansatz to illustrate the behavior due to point 
nodes. These gap functions are illustrated in the insets of FigsJH (a) and (b). Here the 
clean limit ¥ <C -§- is assumed, where T is the carrier scattering rate. 

According to the Doppler shift picture, there are two major differences in the 
angular variation of the thermal conductivity between point nodes and line nodes. First, 
the shape of the k(4>) curves is different. While the oscillation is close to a sinusoidal 
wave for line node (Fig. H](b)), a narrow cusp structure is predicted for the point node at 
T = (Fig. Hta)). Qualitatively, the cusp appears as a result of the small phase space 
available for the quasiparticles induced in the vicinity of point nodes by the applied 
field. For line node the corresponding phase space is greater, and the minimum is not 
as sharp. Second, the amplitude of the oscillation at T = decreases rapidly when 
the H is rotated conically as a function of <fi keeping 9 constant. For point nodes, the 
amplitude of the oscillation of the thermal conductivity at 9 = 45° is much smaller than 
that at 9 = 90°, while they are of almost the same magnitude for line nodes. This can 
be accounted for considering the fact that for 9 = 45° geometry the field H is never 
aligned with the point nodes on the equator. Hence there is always a finite Doppler shift 
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at all the nodes. In contrast, for vertical line nodes, the rotating H at any 9 always 
crosses the line of nodes leading to a greater suppression of the DOS. 

While no microscopic calculations exit at present for this geometry, it is likely that 
the main conclusions of the Doppler shift picture remain valid. As discussed above, 
the salient features of the measurement in the conical experimental geometry are less 
sensitive to the vortex scattering than those measured with the heat current in the 
basal plane. It is likely that the sharp cusp is smeared by finite temperature, but the 
rapid decay of the oscillations as the field is tilted away from the plane must remain 
observable. Nonetheless, more work utilizing microscopic theory is clearly desirable in 
this situation. 

In summary, the variation of the field direction in (9, 0) leads to periodic variations 
in both the thermal conductivity and the heat capacity of nodal superconductors. From 
the periodicity, phase and shape of the angular variation of the thermal conductivity 
and heat capacity, one can extract information on the direction and type of nodes in 
k-space. 

5.3. Experimental 

In the experiments described below the thermal conductivity was measured in a 3 He 
cryostat by using the standard steady-state method, with a heater and two carefully 
calibrated Ru02 thermometers. In all the measurements, and especially in quasi-2D 
superconductors with very anisotropic upper critical field, it is critically important to 
align H in the plane with high accuracy, and have a good control over its rotation. Even 
a slight field-misalignment may produce a large effect on the measured k, influencing the 
conclusions. To achieve this high precision for the orientation of H relative to the crystal 
axes, we used a system with two superconducting magnets generating magnetic fields 
in two mutually orthogonal directions, and a 3 He cryostat set on a mechanical rotating 
stage at the top of a Dewar. By computer-controlling the two superconducting magnets 
and rotating stage, we were able to rotate H with a misalignment of less than 0.02° 
from each axis, which we confirmed by simultaneous measurements of the resistivity. 

Since the thermal conductivity can be measured both under the field 
H=H (sin 9 cos <j), sin # sin 0, cos#) rotated within the basal afr-plane (as a function of 
at 9 = 90), and H rotated as a function of 9 at fixed 0, we were able to detect both 
vertical and horizontal nodal structure. In addition, measuring the thermal conductivity 
with H=H (sin 9 cos 4>, sin sin 0, cos 0) rotated conically as a function of 0, keeping 9 
constant, as shown in FigJU enables us, at least in principle, to distinguish line and point 
nodes. In the following we discuss the experimental results for different compounds. 

One of the recurring aspects of the discussion is the relative importance of electron 
and phonon (or spin-wave) contributions to the net thermal conductivity. While at low 
temperatures the bosons are less efficient than fermions in carrying heat, in systems 
with significant spin fluctuations or low carrier density the bosonic degrees of freedom 
may be dominant over a wide T-H range. Since only electrons carry charge current, 
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we take the point of view that if the measured Wiedemann-Franz ratio of the thermal 
to electrical conductivities, L = n zz p zz /T, just above T c is close to the Lorenz number 
Lq = 2.44x 10~ 8 fiW/K, obtained under assumption of purely electronic k, the electronic 
contribution to the thermal conductivity is dominant. Even though, strictly speaking, 
in the presence of inelastic scattering, L — ► L only as T — > 0, quite generally opening 
of an additional bosonic conduction channel increases L compared to L , and has clear 
experimental signatures. 

In several of the systems we discuss k(T)/T > k(T c )/T c at least in some range 
T < T c and often has peak at T* < T c . In compounds with low carrier density, and 
therefore with thermal transport dominated by bosonic degrees of freedom, this may be 
due to increased mean free path of phonons as the unpaired electron density is decreased. 
In other materials, such increase is due to rapid reduction of the inelastic scattering rate 
below T c (faster than the concomitant reduction in the electron density of states). In 
some systems, both effects combine. Generally, a comprehensive analysis of a large body 
of data on a given compound provides a clue to what mechanism is more important. 
We indicate this for each of the materials analysed below. 

6. Three Dimensional Unconventional superconductors 

6.1. Borocarbide YNi 2 B 2 C 

We start by considering the superconducting gap structure of a non-magnetic 
borocarbide superconductors LnNi 2 B 2 C, Ln=(Y and Lu)[12]. These systems have 
tetragonal crystal symmetry, and the electronic band structure is essentially 3D, see 
Figj5j Early on, these materials were assumed to have an isotropic s-wave gap, similar 
to most compounds where superconductivity is mediated by conventional electron- 
phonon interactions. However, recent experimental studies, such as specific heat [88, 89J, 
thermal conductivity [90], Raman scattering [9~T] . and photoemission spectroscopy [92] 
on YM2B2C or LUM2B2C have reported a large anisotropy in the gap function. Below we 
review the implications of the thermal conductivity measurements for the gap symmetry. 

Figure [6] (a) shows the T-dependence of the c-axis thermal conductivity k zz (the 
heat current q || c) of YNi 2 B 2 C (T c =15.5 K) single crystal with no magnetic field. The 
residual resistivity ratio of this crystal is approximately 47 (the highest crystal quality 
currently achievable). Upon entering the superconducting state, k zz exhibits a small 
kink, as expected for a second order transition. The Wiedemann-Franz ratio at T c , 
L = n zz p zz /T ~ 1.02Lo indicating that the electronic contribution to k is dominant. 
The inset of Figj6] (a) shows the same data below 1 K, where the T-dependence of k zz is 
close to quadratic (rather than cubic, as it would be for dominant phonon contribution). 
Figure[6] (b) depicts the magnetic field dependence of k zz (H\\ [110]) at low temperaures. 
Rapid increase of k zz at low fields is markedly different from that observed in typical s- 
wave materials [54J . This steep increase of the thermal conductivity, along with the \J~H - 
dependence of the heat capacity [881 EH] , strongly suggests that the thermal properties 
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are governed by the delocalized QPs arising from the nodes (or extremely deep minima) 
in the gap. 

Having established the predominant contribution of the extended QPs in the 
thermal transport, we are in the position to address the nodal structure of the gap 
function. As discussed above, in three-dimensional systems, conical rotation of the field 
allows a more direct observation of the nodal structure. Figure [7J displays the angular 
variation of k zz , measured by rotating H=H (sin8 cos 0, sin9 sin 0, cosd) conically, as a 
function of 0, at a constant 9. The measurements were done by rotating after field 
cooling at = —45°. The open circles in FigJTJ show K zz (H,<p) at H—l T which are 
obtained under the field cooling at each angle, and demonstrate excellent agreement 
between the two sets of measurements. A clear fourfold symmetry is observed for the 
0- rotation at 6 = 90° and 60°, so that k zz = k° zz + Here k° zz is 0- independent, and 
fc^j has the fourfold symmetry with respect to 0-rotation. 

As seen in Fig. [7J nff has a narrow cusp at = 0° and 90°. We stress that the 
anisotropies of the Fermi velocity Vp and H C 2, which are inherent to the tetragonal band 
structure of YNi 2 B 2 C, are unlikely to be at the origin of the observed fourfold symmetry. 
The 4-fold 0-dependence of H c2 at 9 = 90° and 45° is nearly perfectly sinusoidal |93j, 
and therefore different from the 0-dependence of k^ z displayed in Fig. [7J According to 
the previous section, the minima of k^ z at = 0° and 90° immediately indicate that the 
nodes are located along [100] and [010] -directions. 

The cusp structure and the ^-dependence of K^f are key features for specifying 
the type of nodes. First, the cusp itself is markedly different from the smooth (almost 
sinusoidal) feature predicted (see previous section, FigH]) and observed (see next section) 
in superconductors with line nodes, such as c/-wave. Second, the amplitude of K,f^ 
decreases rapidly as H is changed from the in-plane 9 = 90° to 45°. Therefore 




Figure 5. Crystal structure of YN^E^C. 
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Figure 6. (a) Temperature dependence of the c-axis thermal conductivity k zz in zero 
field. Inset: Log- log plot of the same data below IK. (b) Field dependence of n zz at low 
temperatures (H\\ [110]). The solid circles represent the data measured by sweeping H 
after zero field cooling, and the open circles represent the data measured under field 
cooling conditions at each temperature. 
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Figure 7. Angular variation of k Z2 (q|| c), measured by rotating H(9,<fi) — 
H (sin 9 cos (j>, sin 9 sin cf>, cos 9) conically as a function of </> at fixed 9 — 90°, 60°, and 
45° (see the inset). The open circles represent the data obtained under the field cooling 
condition at each angle 



direct comparison of the data on both the cusp structure and ^-dependence of k zz 
with FigHJ) strongly favors a model with point nodes, and leads us to conclude that 
the superconducting gap function of YN12B2C has point nodes along [100] and [010]- 
directions. 

For a gap with point nodes the T-dependence of the thermodynamic quantities at 
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low temperature depends on whether the gap increases linearly or quadratically with 
the distance from the nodal point. The gap function we used predicts a quadratic T- 
dependence of the thermal conductivity, which is consistent with the data in the inset of 
Fig. W(a). 

More recent measurements of the angular variation of the heat capacity also report 
the four fold oscillations consistent with the present experiments [70] . Very recent STS 
measurements in the vortex state of YNi 2 B 2 C have demonstrated the presence of the 
extended QPs in the [010] direction [91] . Thus the nodal structure is confirmed by the 
several different techniques. 

Are these real nodes or simply deep minima? Experimentally, a clear fourfold 
pattern is seen at T=0.27 K and H = 1T~ 0.1H c2 . This suggests that the typical 
Doppler energy, E av ~ A a/ H/H c2 , of the nodal quasiparticles far exceeds both T, and 
the minimal gap A min . Here we estimate A ~ hvp/n^Q. This leads to the anisotropy 
ratio A min /Ao -C 0.3. A more stringent constraint may be deduced from the power law 
temperature dependence of the thermal conductivity down to this temperature in zero 
field. Estimating, A ~ 28K from the value of T c , we find A min /A < 0.01. 

While this value is small, the origin of the true nodes is topological, and hence 
the important question is whether the gap function changes its sign on the Fermi 
surface. To answer it, we examined the impurity effect on the gap anisotropy. In an 
anisotropic s-wave superconductor, with accidental gap minima or zeroes, introduction 
of non-magnetic impurities affects T c only moderately, and rapidly makes the gap more 
isotropic thereby reducing the DOS at the Fermi surface by removing the node. On 
the other hand, if the gap changes sign and its average over the Fermi surface vanishes, 
doping with impurities suppresses T c more severely, and induces a finite DOS at energies 
smaller than the scattering rate, 7. In the latter case the oscillations of k(H) persist in 
the regime where the Doppler energy 7 < E av . 

Figure [8] shows the absence of the angular variation of the thermal conductivity 
k zz (6 = 90°, 0) in the Pt-substituted compound Y^^Pi^C with x=0.05. We 
estimate that E av ~ 0.3Ao- On the other hand, the transition temperature changes 
little with Pt-doping [881 [95]. Then the disappearance of the angular variation in k 
indicates the opening of the gap, and the destruction of nodal regions by impurity 
scattering. This is consistent with the heat capacity measurements, which reports a 
transition from \/H behavior of 'y(H) at x=0 to linear in H behavior for x=Q.2 [88] . 

It has been pointed out that the cusp structure in the angular variation of the 
thermal conductivity can appear as a result of the nesting property of the Fermi surface 
[80J. However, the disappearance of the angular variation in Y(Ni 1 _ a ,Pt x )2C with 
£=0.05, indicates that this scenario is unlikely. Moreover the cusp structure appears 
even in ErN^E^C, in which the nesting part of the Fermi surface disappears due to the 
spin-density- wave transition [96] . 

Therefore, comparison of the experiment with existing theories yields the gap 
structure with point nodes. While these may be accidental, an alternative view is 
that a strong Coulomb repulsion is an essential ingredient of the models required for 
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Figure 8. Angular variation of k zz for YM2B2C (open circle) and Y(Nio.95Pto.o5)2B2C 
(solid circles), measured by rotating H as a function of cj> within the afe-plane. Angular 
variation disappears in Y(Nio.95Pto. 05)2820. 



the borocarbides. Recently, it has been shown that the s-wave superconductivity with 
deep gap minimum appears when the electron phonon coupling coexists with the AF 
fluctuation [97J. We note that the topology of the nodal regions plays an important 
role in determining the superconducting properties, such as the vortex lattice structure, 
reversible magnetization, upper critical field etc. For instance, the extended QPs 
appear to be very important for the vortex triangular-square lattice phase transition 

6.2. Heavy Fermion UPCI2AI3 

UPd2Al 3 has aroused great interest among heavy fermion (HF) superconductors because 
of its unique properties. In UPd2Al 3 , superconductivity with heavy mass occurs 
at T c =2.0 K after antiferromagnetic (AF) ordering with atomic size local moments 
(p =0.85/is) sets in at T/v=14.3 K [100J. Below T c , superconductivity coexists with 
magnetic ordering. The ordered moments are coupled ferromagnetically in the basal 
hexagonal afe-plane and line up along the a-axis (Fig. 9). These ferromagnetic sheets are 
stacked antiferromagnetically along the c-axis with the wave vector Q o =(0,0,ir/c), where 
c is the c-axis lattice constant |101j (For the structure of the hexagonal basal plane, see 
the inset of FigJTQl The presence of large local moments is in contrast to other HF 
superconductors, in which static magnetic moments are either absent or very small at 
the onset of superconductivity [6j. Since both superconductivity and AF ordering in 
UPd2Al 3 involve the Uranium 5/ electrons, this system is an example of the dual nature, 
partly localized and partly itinerant, of strongly correlated electrons [6l 11021 [103] . 

In the superconducting state of UPd2Al3, two noticeable features have been 
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reported. The first is the "strong coupling anomalies" observed in the tunnel junctions 
UPd2Ai 3 -A10. 2 ;-Pb 1 104] and inelastic neutron scattering, which were attributed to 
the strong interaction between the heavy fermion quasiparticles and AF spin-wave 
excitations [1031 1105[ I106[ I107j . The second feature is the appearance of a "resonance 
peak" in the inelastic neutron scattering in the vicinity of Qo well below T c . Similar 
peak in the cuprates was interpreted as the result of the feedback effect of the opening of 
the superconducting gap on the electron-hole damping of the spin fluctuations, and was 
shown to be unique to superconductors where the gap changes sign under the translation 
by the AFM wave vector [108J. The analogous effect in UPd2Al 3 , which has static 
AFM order, rather than fluctuations, was investigated by Bernhoeft et al. [105l I106j . 
and strongly suggests that the gap in this material changes sign under translation k— > 
k + Qo- The NMR Knight shift measurements indicate the spin-singlet pairing, and 
the spin-lattice relaxation rate does not show the coherence peak at T c , and decreases 
as Tjf oc T 3 , indicating the presence of line nodes |109j . These results provide rigorous 
constraints on the shape of the gap. 

According to band calculations and de Haas-van Alphen measurements in the AF 
phase, the largest Fermi sheet with heavy electron mass and the strongest 5f-admixture 
has the shape of a corrugated cylinder with a hexagonal in-plane anisotropy [102J. Below 
we assume that this cylindrical Fermi sheet with heavy mass is responsible for the 
superconductivity, and carry out the analysis within this model. 

We measured the thermal conductivity along the c-axis of the hexagonal crystal 
structure, k zz (heat current q || c) and along the 5-axis n yy (q \\ b) in high quality single 
crystals of UPd2Al 3 with T c = 2.0 K (The residual resistivity ratio was 55 along the 
6-axis and 40 along the c-axis). Figure [TOl depicts the temperature dependence of n yy 
and k zz in zero field. Since spin- wave spectrum has a finite gap of ~ 1.5 meV at the zone 
center, its contribution appears to be negligible below T c [109[ 1110] . The Wiedemann- 
Franz ratio L = at T c is 0.95Lq for k zz and is 1.16-Lq for K yy . These results indicate 
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Figure 9. Crystal structure of UPCI2AI3. AF ordering with atomic size local moments 
(/j =0.85^) sets in at Tjv=14.3 K. The ordered moments are coupled ferromagnetically 
in the basal hexagonal afr-plane and line up along the a-axis. These ferromagnetic 
sheets are stacked antiferromagnetically along the c-axis. 
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Figure 10. Temperature dependence of the thermal conductivity along the &-axis K yy 
(q|| b) and c axis k zz (q|| c)in zero field. Inset: structure of the hexagonal basal plane 
of UPd 2 Al 3 with the alignment of the a-axis (100) and 6-axis (-1,2,0), as used in the 
text. 
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Figure 11. Field dependence of the 6-axis thermal conductivity K yy normalized to 
the normal state value just above the upper critical field, ft" , for H\\ a and H\\ c. 



that the electron contribution is dominant below T c . 

Figure [TT] shows the if- dependence of K yy for H\\ a and JEf|| c below T c . For both 
field directions, K yy grows with H beyond an initial decrease at low fields. For H || c, 
K yy increases almost linearly with H, K yy oc H, at 0.36 K. The minimum in k(H) is 
much less pronounced at lower temperatures. As H approaches H c2 || a or || b, K yy 
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Figure 12. Angular variation of the c-axis thermal conductivity n zz (H,cf)) normalized 
by the normal state value, k" z , at several fields at 0.4 K, below i/ C 2(=3.5 T). H was 
rotated within the basal plane as a function of <f> (see the inset). A distinct six- fold 
oscillation is observed above 0.5 T, while oscillation is absent at 0.5 and 0.3 T. 



shows a steep increase and attains its normal state value. This low to intermediate 
if -dependence of the thermal conductivity in the superconducting state is markedly 
different from that observed in ordinary s-wave superconductors. At high temperatures 
and low fields, where the condition a/ H/ H c2 < T/T c is satisfied, the thermally excited 
quasiparticles dominate over the Doppler shifted quasiparticles. It has been shown 
that in this regime, while the Doppler shift enhances the DOS, it also leads to a 
concomitant reduction in both the impurity scattering time and Andreev scattering 
time off the vortices [32J EH EU [65j 1111] . When this lifetime suppression exceeds the 
enhancement in N(E), which may happen at intermediate temperatures and low fields, 
the nonmonotonic field dependence of the thermal conductivity is found. As in other 
superconductors with nodes, the region of the initial decrease of the thermal conductivity 
shrinks at low T. Thus the if -dependence of K yy in UPd2Al3, initial decrease at low 
field at high temperatures and linear behavior K yy oc H at low temperatures, are in 
qualitative agreement with the existence of line nodes in A (JET) [82j 1112] , 

We first test whether there exist vertical line nodes perpendicular to the basal 
plane. Figure [T2l shows K zz (H,<j)) as a function of <fi at 0.4 K, measured by rotating H 
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Figure 13. Angular variation of the &-axis thermal conductivity K yy (H,9) normalized 

by the normal state value k™ at several fields at 0.4 K. H was rotated within the 

J yy 

ac-plane perpendicular to the basal plane (see the inset). 

within the basal plane (9 = 90°). Above 0.5 T, a distinct six-fold oscillation is observed 
in K zz (H,(f)), reflecting the hexagonal symmetry of the crystal. Six- fold oscillation is 
observable even above H C 2. On the other hand, no discernible six- fold oscillation was 
observed below 0.5 T within our experimental resolution. We infer that the AF magnetic 
domain structure and anisotropy of H c2 within the plane are responsible for the six-fold 
symmetry and the nodal structure is not related to the oscillation. According to the 
neutron diffraction experiments, the magnetic domain structure changes at Hp ~0.6 T 
well below T/y. Below Hp, the ordered moments point to the a-axis, forming domains, 
and the spin structure is not affected by the ff-rotation in the basal plane. On the 
other hand, above Hjj, the H -rotation causes domain reorientation. Then the magnetic 
domain structure changes with sixfold symmetry with //-rotation. Thus the sixfold 
symmetry observed in k zz above 0.5 T is most likely to be due to the magnetic domain 
structure. This indicates that there are no nodes located perpendicular to the basal plane, 
i.e. gap function in the basal plane is isotropic. 

We next test for the existence of the horizontal line nodes parallel to the basal plane. 
Figure [131 displays the angular variation of K yy (H,9) for rotating H as a function of 9 
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within the ac-plane at 0.4 K. A distinct oscillation with two-fold symmetry was observed 
in the superconducting state. In contrast to k zz , no discernible twofold oscillation was 
found in the normal state above H c2 . We decompose K yy (H,9) as 



where n yy is a independent term and K yy = C yy cos 29 is a term with the twofold 
symmetry with respect to ^-rotation. The field dependence of C yy at 0.4 K is shown 
in Figure [141 (a). For comparison, the independence of K yy for ff|| c at T=0.36 K is 
plotted in Figfl4l (b) . There are three regions denoted (I), (II) and (III), below H c2 . In 
the vicinity of H c2 ((Ill)-region), where n yy increases steeply with H, the sign of C yy 
is negative and the amplitude \C yy \/n yy is of the order of 10%. Here n yy is K yy in the 
normal state just above H c2 . With decreasing H, C yy changes sign at about 2.3 T and 
becomes positive in the region where K yy for H || c shows a linear if -dependence ((H)- 
region). Below about 0.25 T, where the second sign change takes place, n yy decreases 
with H ((I)-region). In this region, C yy is, once again, negative. 

We address the origin of the observed two-fold oscillation. The disappearance of 
the oscillation above H c2 , together with the fact that there is only one magnetic phase in 
this configuration [101J, completely rule out the possibility that the origin is due to the 
magnetic domain structure. There are two possible origins for the oscillation; the nodal 
structure and the anisotropy of the Fermi velocity and H c2 . Obviously, as discussed 
previously, a large two-fold oscillation with negative sign observed in the (Ill)-region 
arises from the anisotropics of the Fermi velocity and H c2 . This immediately indicates 
that the two-fold symmetry with positive sign in the (Il)-region originates not from these 
anisotropics but from the quasiparticle structure associated with the nodal gap function. 
In addition, the amplitude of C yy /n yy in the (Il)-region is a few percent, which is 
quantitatively consistent with the prediction based on the Doppler shifted DOS. We 
also note that the second sign change at low fields in the (I)-region is compatible with 
the nodal structure. In this region, as discussed previously, the if -dependence of the 
thermal conductivity is governed by the suppression of the quasiparticle scattering rate. 

As discussed in Ref . [Ml [58j EH ED EU 1111] , the anisotropic carrier scattering time 
associated with the nodal structure also gives rise to the variation of K yy (H,9) as a 
function of 9. In this case the sign of the oscillation is opposite to that arising from the 
Doppler shifted DOS in the (Il)-region. These considerations lead us to conclude that 
UPd2Al 3 has horizontal nodes. In addition, the fact that there is a single maximum 
structure in the angular variation of n yy (H,9) indicates that horizontal line nodes are 
located at positions where the condition p|| ab in the Brillouin zone is satisfied. Thus 
the allowed positions of the horizontal nodes are restricted at the bottleneck and AF zone 
boundary [621 182] . 



For comparison, the angular variations of K yy and k zz at low fields are shown 
in FigJT5l While the amplitude of the two-fold oscillation C yy /n yy is 3%, which is 
quantitatively consistent with the Doppler shifted DOS, the amplitude of the six-fold 
oscillation C® z /k zz is less than 0.2%, which is more than 10 times smaller than the 
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Figure 14. (a) Field dependence of the amplitude of the two-fold symmetry C yy 
normalized by the normal state thermal conductivity K yy at 0.4 K. (b) Field dependence 
of the 6-axis thermal conductivity n yy for ff|| c at 0.36 K. The sign of C yy is negative 
in the (Ill)-region just below H C 2 and in the (I)-region where K yy decreases with H . 
On the other hand, the sign of C yy is positive in the (Il)-region. 



amplitude expected from the Doppler shifted DOS in the presence of nodes. Combining 
the results, we arrive at the conclusion that the gap function is isotropic in the basal 
plane and has horizontal node. 

To discuss the position of the horizontal line node in UPd2Al3, we consider a 
"magnetic" Brillouin zone in a cylindrical Fermi surface, as shown in Fig{16] (a). The 
density of states is anisotropic under the rotation of the field in the ac-plane, by varying 
the angle 9 in the inset. To illustrate the difference between the line nodes at high 
symmetry positions in the magnetic Brillouin zone, and away from those, we consider 
here four model gap functions 

(i) type-I : A horizontal node located at the bottleneck; A(k) oc sin k z c. 

(ii) type-II : A horizontal node located at the zone boundary; A(k) oc cosfc^c. 

(iii) type-Ill: A hybrid of type-I and -II. Two horizontal nodes located at the bottleneck 
and the zone boundary; A(k) oc sin2/c 2 c. 

(iv) type-IV : Two horizontal nodes located at positions shifted off the bottleneck in 
the Brillouin zone; A(k) oc cos2k z c. 

The expected angular variation of the Doppler shifted DOS is a function of the relative 
angle between H and p for these gap functions are shown schematically in Fig{16] (b). 
The twofold oscillations with the same phase are expected for type-I, -II , and -III gap 
functions, in which the horizontal nodes are located at the position where p|| afr-plane; 
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Figure 15. Angular variation of the 6-axis thermal conductivity k^/k" with rotating 
H within the ac-plane (solid circles) and of the c-axis thermal conductivity k^/k™ z 
with rotating H within the basal a6-plane (open circles) at T=0.4 K and at H=0.5 T. 
The amplitude of the six-fold oscillation in K z f/K^ z is less than 0.2% if it exists. 
Inset: Schematic figure of the gap function of UPd2Al3 determined by angle resolved 
magnetothermal transport measurements. The thick solid lines indicate horizontal 
nodes located at the AF zone boundaries. 



one cannot distinguish these three gap functions when the Fermi surface has an open 
orbit along the c-axis. For type-IV, one expects an oscillation with a double minimum 
structure as a function of 6. 

Thus, the order parameters allowed, by thermal conductivity measurements, are, 

(i) A(k) = A sinA; 2 c, 

(ii) A(k) = Aosin2A; z c and 

(iii) A(k) = Aocos/c^c. 

which are shown in the type-I, -II and -III gap structures in Fig lTBlf a). Generally the 
first and the second represent spin triplet gap functions, and only the third, which is a 
spin singlet, remains a viable possibility. Note, however, that the thermal conductivity 
is only sensitive to the amplitude of the gap. Among the possible order parameters 
for the D 6 h symmetry group is that transforming according to T 5 representation, 
with a basis function k z (k x + ik y ) pQ. In that case the gap function may vary as 
) sin k z c } and the gap amplitude, over a quasi-two dimensional Fermi surface, 
may have only weak modulations apart from the horizontal line of nodes. Such an order 
parameter breaks the time reversal symmetry, and we are not aware of any evidence 
in support of that in UPd2Al3; however, targeted search for a time-reversal symmetry 
broken state in this system has not been performed. Both the first and the third gap 
functions are compatible with the constraint implied by the neutron resonance peak, 
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Figure 16. (a) Schematic figure of the gap structure with horizontal line node in 
the magnetic Brillouin zone. Line nodes are located at the bottleneck (type-I) and 
at the zone boundary (type-II). Two line nodes are located at the bottleneck and the 
zone boundary (type-Ill) and at positions shifted off the bottleneck (type-IV) . (b) 
Oscillations of the DOS for H rotating in the ac-plane for various gap functions. Two- 
fold oscillation with the same sign are expected for type-I, -II, and -III. On the other 
hand, for type-IV, an oscillation with a double minimum is expected. 



A(k) = — A(k+Q ). These considerations lead us to conclude that the gap function of 
UPd 2 Al 3 is most likely to be A(k) = A cos k z c, shown in the inset of Figfl~5l although 
we cannot exclude the possibility of the state with broken time-reveral symmetry on the 
basis of out measurements. 

For the cos k z c pairing, the horizontal node located at the AF zone boundary 
indicates that pair partners cannot reside in the same basal plane. The interlayer 
pairing appears to indicate that strong dispersion of the magnetic excitation along k z 
causes the pairing, as suggested in the magnetic exciton mediated superconductivity 
model P, 11131 1103] . The isotropic gap function in the basal plane implies that the 
pairing interaction in the neighboring planes strongly dominates over the interaction 
in the same plane. Although the pairing interaction inferred from the determined gap 
function should be further scrutinized, the recent results imply that the interlayer pairing 
interaction associated with the AF interaction is most likely to be the origin of the 
unconventional superconductivity in UPd 2 Al 3 . 

6.3. skutterudite PrOs^Sb^ 

Recently discovered heavy fermion superconductor PrOs4Sbi2 (T c =1.82 K) with filled 
skutterudite structure (Fig. 17) is relatively unique as the /-electrons have a non- 
magnetic ground state, determined by the crystalline electric field (most likely singlet T\ 
state) [114[ 1115] . The HF behavior (m* ~ 50m e , m e is the free electron mass) is likely 
due to the interaction of the electric quadrupole moments of Pr 3+ (rather than local 
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magnetic moments as in the other HF superconductors) with the conduction electrons. 
Therefore the relation between the superconductivity and the orbital (quadrupole) 
fluctuations of /-electron state excited great interest: PrOs4Sbi2 has been proposed as a 
candidate for the first superconductor with pairing mediated neither by electron-phonon 
nor magnetic interactions, but by quadrupolar fluctuations. Even if these fluctuations do 
not provide the pairing glue by themselves, but only in conjunction with phonons, they 
have the potential for influencing the symmetry of the superconducting order parameter, 
which makes it of the utmost importance to determine the symmetry of the SC gap. 

The unconventional superconductivity in PrOs 4 Sbi 2 has been suggested by several 
experiments. NQR measurements showed the absence of Hebel-Slichter peak [116J. In 
addition, the Knight shift does not change below T c , implying that the gap function 
has odd parity [22J. Moreover, /zSR experiments report the appearance of the static 
spontaneous magnetic field below T c , which can be interpreted as the spontaneous 
breaking of the time reversal symmetry [33J. The penetration depth and NMR T{~ 
measurements indicate the presence of point nodes [117] . 

Figure [18] shows the T-dependence of the c-axis thermal conductivity k zz (the heat 
current q || c) divided by T both at zero field and above H c2 (~ 2.2 T at T— K) of 
PrOs 4 Sbi2 single crystal (Jjl). In this temperature region, the electronic contribution to 
k zz dominates the phonon contribution. The inset of Fig{18] shows the field dependence 
of k zz of sample Jj2 at very low temperature. k zz increases very steeply even at very low 
field (H < 0.1T). When contrasted with the exponentially slow increase of the thermal 
conductivity with field observed in s-wave superconductors at H -C H c2 |118j . this is a 
strong indication that the thermal transport is governed by the delocalized QPs arising 
from the gap nodes. Above 0.1 T k zz increases gradually, then shows a steep increase 
above 0.5 T up to H c2 . 

Figures [191 (a) and (b) display the angular variation of K zz (H,<p) in H rotated 
within the a6-plane (9 = 90°) at T=0.52 K [68]. The measurements have been done in 
rotating H after field cooling at <fi = —90°. The open circles show k zz (H, 0) at if =1.2 T 




Figure 17. Crystal structure of PrOs4Sbi4 
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Figure 18. Temperature dependence of the c-axis thermal conductivity k zz (the heat 
current q || c) divided by T at zero field (solid circles) and at 2.5 T (open circles) above 
H c2 (~ 2.2 T at T= K) of PrOs 4 Sbi 2 single crystal (jjl). Inset: Field dependence of 
Kzz/T of sample j}2 at very low temperature. 
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Figure 19. (a)(b) Angular variation of n zz (q|| c) in H rotating within the a6-planc 
as a function of 4> at 0.52 K above and below H c2 {~ 2.0 T). 



and 0.5 T which are obtained under field cooling at each angle. Above H C 2 (~2.0 T at 
0.5 K) K zz (H,(f>) is essentially independent of 0. A clear fourfold variation is observed 
just below H c2 down to H ~ 0.8 T. However further reduction of H below 0.8 T causes 
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Figure 20. The amplitude of twofold (open circles) and fourfold (filled circles) 
symmetries, |C 2 ^| and \C^\, respectively, plotted as a function of H/H c2 at T=0.52 K. 
At H*, the crossover from twofold to fourfold symmetry takes place. 



a rapid decrease of the amplitude of the fourfold term, and its disappearance below 
0.7 T. At the same time, the twofold component grows rapidly. This surprising behavior 
suggests a change in the gap symmetry as a function of field in the superconducting state. 

Figure [201 shows the if -dependence of the amplitudes of the twofold and the fourfold 
terms, which are obtained by decomposing K zz (H,(f)) as 

k kk {H,4) = i& + i& + i£, (14) 

where n° zz is a 0-independent term, k 2 £ = C 2 cos 20, and = cos 40 are the terms 
with twofold and fourfold symmetry with respect to 0-rotation. It is clear that the 
transition from the fourfold to twofold symmetry in 0-rotation is sharp, and occurs in a 
narrow field range at H/H C 2 ~ 0.4, deep inside the SC phase. Both symmetries coexist 
in a narrow field range. If the minima of the thermal conductivity are associated with 
the direction of the field along the nodes, the reduced K zz (H,<p) at = ±90° and 0° 
in the high field phase, and at = ±90° for the low field phase, respectively, lead to 
the conclusion that the nodes are along the [100]- and [010]-directions in the high field 
phase, while they are located only along the [010]-direction in the low field phase. 

Having established the presence of nodes, the next question is their classification. 
As discussed in § IV A, the angular variation of 

f^ZZ Colli distinguish between the point 
and line nodes, by rotating H conically around the c-axis with a tilted angle from the 
afr-plane. Although we do not show it here, the amplitude at 9 = 45° and 30° is smaller 
than that at 9 = 90° [HE]- Similar results were obtained for the twofold symmetry. 

What is the nodal structure inferred from the present results? The 0-rotation of 
the field can only provide the information of the nodes away from the [001] direction. 
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Figure 21. The phase diagram of the superconducting gap symmetry determined by 
the present experiments. The hlled circles represent the magnetic field H* at which the 
transition from fourfold to twofold symmetry takes place. The open circles represent 
H C 2- The area of the gap function with fourfold symmetry is shown as A-phase and 
the area of the gap function with twofold symmetry is shown as B-phase. 



We shall therefore appeal to the group theoretical consideration for the discussion of the 
nodal structure. It is unlikely that the SC gap function has only four point nodes in the 
cubic Th crystal symmetry: this is independent of the spin singlet or triplet symmetry. 
Hence the likely scenario is that the gap function at high field phase has six point nodes. 
The low field phase is likely to have two nodes, although on the basis of the experimental 
observations we cannot exclude the 4-node structure (along [001] and [010]) 

The H — T phase diagram of the SC symmetry determined by the present 
experiments is displayed in Fig. [2TJ The filled circles represent the magnetic field 
H* at which the transition from fourfold to twofold symmetry takes place. The H*- 
line which separates two SC phases (high field A-phase and low field 5-phase) lies 
deep inside the SC state. We note that recent flux flow resistivity measurements also 
reported an anomaly at H*. The only example of a superconductor with multiple phases 
of different gap symmetry so far has been UPt3 |119j . In that case the degenerate 
transition temperatures for the two orders at zero field can be split by, for example, 
applying pressure. Therefore it is important to determine a) whether the two phases 
manifested in the thermal conductivity measurements have the same T c in zero field; 
b) whether the transition can be split by influencing the system by an experimental 
handle other than the field. It seems logical that, if the gap structure suggested here is 
indeed realized, application of the uniaxial pressure along the [100] direction should lift 
the symmetry of the gap and favor one of the two phases. 

Recently, small angle neutron scattering experiments reported the hexagonal flux 
line lattice, which is distorted with a twofold symmetry |120j . It has been pointed out 
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that the distortion originates from the nodal gap structure, which provides a strong 
support of the present angular variation of K zz (H,<ft). We also note that the possible 
existence of the third phase was predicted by the magnetization and penetration depth 
measurements [1171 1121 j . 

Thus PrOs4Sbi2 has several unique features. In almost all superconducting (SC) 
materials known to date, once the energy gap in the spectrum of electrons opens at 
the SC transition, only its overall amplitude, and not the shape and symmetry around 
the Fermi surface, changes in the SC phase [122j . In contrast, PrOs4Sbi2 seems to 
have several superconducting phases with different symmetries. Many heavy fermion 
superconductors have line nodes in the gap functions (with possible additional point 
nodes). The suggestion that PrOs 4 Sb 12 is the first heavy fermion superconductor, 
in which only the point nodes are identified may be the key for understanding the 
superconductivity mechanism due to quadrupolar interaction. 

7. Quasi two dimensional superconductors 

We now discuss the implications of the thermal conductivity measurements for the nodal 
structure of three superconductors, CeCoIn 5 , k-(BEDT-TTF) 2 Cu(NCS)2, and Sr 2 Ru04. 
These materials look very different at first sight but reveal several similar features. First, 
all three have strong electron-electron correlations. Second, they all have quasi-two 
dimensional electronic structure, as confirmed by the band structure calculation and by 
the dHvA measurements. This is also supported by a relatively large anisotropy of the 
upper critical field between inequivalent crystalline directions. Third, the power laws in 
the temperature dependence of thermodynamic quantities in the superconducting state 
is consistent with the presence of line nodes in the superconducting gap. The position 
of the line nodes is still an open question in many of these materials, and is the focus 
of our analysis here. 

Unfortunately, in all these compounds, the out-of-plane thermal conductivity 
is difficult to measure due to thin plate-like single- crystal samples. We therefore 
measured the in-plane thermal conductivity with magnetic field rotated within the 
same conducting plane. In this geometry, the dominant signal is twofold symmetric, 
and simply depends on the angle between the thermal current, q, and H due to the 
difference in the transport along and normal to the vortices. This twofold oscillation is 
not directly related to the nodal structure, and the challenge for the interpretation is to 
separate it from the features due to the nodes. 

7.1. CeCoIn 5 

The family of the heavy fermion superconductors CeTIn 5 (T=Rh, Ir, and Co) 
was discovered in 2001 [123J. Both Celrlns and CeCohis are ambient pressure 
superconductors, with transition temperatures of 0.4 K and 2.3 K, respectively. CeRhIn 5 
is an antiferromagnet, but shows superconductivity under moderate pressure. The 
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crystal structure of CeTIn 5 is tetragonal consisting of the conducting Celn 3 layers 
separated by less conducting TIn 2 layers (Fig. 22). In the normal state CeCoIn 5 
exhibits non-Fermi-liquid behavior, likely related to the strong AFM fluctuations; 
moreover, there is evidence that superconductivity appears in the neighborhood of a 
quantum critical point (QCP), possibly of AFM origin pSl EEl [1261 HZZj. For that 
reason CeCoIn 5 is an excellent candidate to study the relationship between QCP and 
unconventional superconductivity. The temperature dependence of the heat capacity 
and thermal conductivity [128J, thermal Hall conductivity [129J, NMR relaxation rate 
|130] all indicate the presences of line nodes. The penetration depth measurements 
|131[ 11321 I133j in general support this conclusion although the low-T exponents are 
anomalous. Recent Andreev reflection measurements indicate the sign change of the 
superconducting order parameter [1341 11351 [136] , see, however |137j for comments on 
Ref. [135j . The Knight shift measurement of CeCoIns indicates the even-spin parity 



There are indications that the upper critical field in CeCoIns is paramagetically 
(Pauli) limited. At low T the phase transition at H c2 is first order, as revealed by a 
step in the if-dependence of the thermal conductivity, magnetization, and specific heat 
[6T | I138[ 1139] . Moreover, measurements of heat capacity [140[ 1141] . ultrasound [142] . and 
NMR |143j revealed a new superconducting phase at low temperatures in the vicinity of 
H c2 at low temperatures (H || ab). This new phase was conjectured to be the spatially 
inhomogeneous superconducting state (Fulde-Ferrell-Larkin-Ovchinnikov state), which 
was predicted 4 decades ago, not previously observed. 

The inset of Fig ]237 a) shows the T-dependence of k and p. Upon entering the 
superconducting state, k exhibits a sharp kink and rises to the maximum value at 
T ~ 1.7 K. The Wiedemann- Franz ratio L = ^p ~ 1.02Lo at T c is very close to the 
Lorenz number Lq = 2.44 x 10 -8 QW/K, indicating that the electronic contribution 
is dominant. Therefore the enhancement of k below T c is due to the suppression of 
the inelastic scattering rate, similar to the high-T c cuprates. Figures [23] (a) and (b) 
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Figure 22. Crystal structure of CeCoIn 5 . 
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Figure 23. Field dependence of the a-axis thermal conductivity (q||[100]) for (a) 
H || [001] and (b) if|| [010] below T c . Inset of (a) : k and p in zero field. Inset of (b) : 
Field dependence of k near H C 2 at 0.45 K (o) and 0.34 K (•). The thermal conductivity 
jumps at H C 2, indicating the first order phase transition. 



depict iJ-dependence of k for ff|| ab (H C 2 — 11 T) and HI. ab (H C 2 — 5 T) below 
T c , respectively. At all temperatures, k decreases with H, and the if- dependence is 
less pronounced at lower T in both configurations. For HA. ab, k has a discontinuous 
jump to the normal state value at H c2 below 1.0 K (see also the inset of Fig j237 b)). 
indicating a first-order phase transition. The data in Figsj23] (a) and (b), in which the 
if-dependence of k is more gradual with decreasing T, are consistent with the picture 
where the scattering of the field-induced QPs is the main origin for the if-dependence 
of K. 

Given the complexity and richness of behavior of this superconductor, it is natural 
to ask whether one can draw conclusions about the symmetry of the gap from the 
analysis of the thermal conductivity rooted essentially in a BCS-like theory. We 
believe that the answer is affirmative. In the following it is important to note that 
all the measurements are done far below the upper critical field, i.e. far away from 
the possible competing states, first order transition, and quantum critical behavior. 
The thermal conductivity is measured at temperatures far below that of the inelastic 
scattering-induced peak in k(T). Consequently, we believe that, similar to the high-T c 
superconductors, the BCS-like model gives the semi-quantitatively correct results in this 
regime. 

Figure [Ml displays n(H,<f)) as a function of <fi =(q, H) at T=0.45 K of CeCoIns 
|61j . The solid circles in FigEMl show K>(H,<f)) at H=l T which are obtained under the 
field cooling at each angle. In all data, as shown by the solid lines in Fig. 2, n(H, <fi) can 
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Figure 24. The a-axis thermal conductivity k(H, <j>) (q|| a) in H rotated within 
the afe-plane as a function of <j> for CeCoIn 5 . 4> is the angle between q and H 
(see the inset). The solid lines represent the result of the fitting by the function 
K,(H,<p) — C + C 2< /,cos2(j) + C40cos40, where C , C 2< p and C±<p are constants. The 
solid circles represent K(H,(f>) at H—l T which are obtained under the field cooling 
condition at every angle. 

be decomposed into three terms, 

K{4>) = K + K 24> + K 44n (15) 

where k is a ^-independent term, and k 2 </> = C 2( p cos 20 and = cos 40 are terms 
with 2- and 4-fold symmetry with respect to the in-plane rotation, respectively. 

Figures [25] (a)- (d) display (normalized by the normal state value K n ). It is clear 
that k 4 exhibits a maximum at if||[110] and [110] at all temperatures. Figure 4 and 
the inset show the T- and H- dependences of |C4^,|/«; n . Below T c the amplitude of 
increases gradually and shows a steep increase below 1 K with decreasing T. At low 
temperatures, {C^l/nn becomes greater than 2%. 

We note that the anisotropy of H c2 (H^ 100 ^ ~ l.OSH^ 1 ^) is too small to explain the 
large amplitude of {C^J K n > 2% at H <^ H c2 . Further, and more importantly, the sign 
of the observed fourfold symmetry is opposite to the one expected from the anisotropy of 
H c2 . The observed 4-fold symmetry above T c is extremely small; \C^\j K n <0.2 %. Thus 
the anisotropies arising from H c2 and the band structure are incompatible with the data. 
Note also that, even if there is a fraction of electrons that remain uncondensed at low 
T, as was recently suggested |144j , they can only influence the twofold (via the orbital 
magnetoresistance), rather than the fourfold symmetry. These considerations lead us 
to conclude that the 4-fold symmetry with large amplitude well below T c originates from 
the QP structure. 

We now address the sign of the 4-fold symmetry. For nodal lines perpendicular 
to the layers, two effects compete in determining k^. The first is the DOS oscillation 
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Figure 25. (a)-(d) The 4-fold symmetry normalized by the normal state value k„ 
at several temperatures. 



under the rotation of H within the a6-plane. The second is the quasiparticle scattering 
off the vortex lattice, which has the same symmetry as the gap function [58], [59[ 1501 [53] . 
As discussed above, it is likely that at T « 0.25T C and H < 0.25if C 2 the second effect is 
dominant [23 E3] • In this case, k attains the maximum value when H is along a node and 
has a minimum when H is directed towards the antinodal directions [1H1 [53], [53 159| 164"]. 
The amplitude of the four-fold symmetry \C^\ in quasi-2D d-wave superconductors is 
roughly estimated to be a few percent of n n [161 [53], which is of the same order as the 
experimental results. 

It is interesting to compare our results on CeCoIn 5 with the corresponding results 
on YBa2Cu307_,5, in which the 4-fold symmetry has been reported in the regime where 
the Andreev scattering dominates. In YBa2Cus07_5 with d x 2_ y 2 symmetry, has 
maxima at Jf||[110] and [110] J58J (59[ [M], in accord with our CeCoIn 5 data. Thus 
the sign of the present fourfold symmetry indicates the superconducting gap with nodes 
located along the (±tt,±tt)- directions, similar to the high-T c cuprates; CeCoIn 5 most 
likely belongs to the d x 2_ y 2 symmetry. 

It is worth commenting on the gap symmetry of CeCoIns determined from the 
other techniques. Small angle neutron scattering experiments have reported the square 
lattice of flux lines |145j . whose orientation relative to the crystal lattice is consistent 
with the expectation for the d x 2_ y 2 -symmetry. Recent point contact spectroscopy, 
which measured the Andreev reflection at the normal metal/CeCons interface with two 
crystallographis orientations, (001) and (110), have also concluded that the symmetry is 
d x 2_ y 2 [136] . In contrast to these results, the angular variation of the heat capacity in H 
rotated within the a6-plane originally was interpreted as evidence for the d xy -symmetTj 
[69] . However, recent theoretical analysis suggested that, when the redistribution of 
the spectral density due to vortex scattering is accounted for, the specific heat is also 
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Figure 26. Crystal structure of k-(BEDT-TTF) 2 Cu(NCS) 2 

consistent with the d x 2_ y 2 gap [53]. 

Therefore most of the measurements suggests that this material has d x i_ y i gap 
symmetry, likely implying that the antiferromagnetic fluctuations are important for 
superconductivity. This observation qualitatively agrees with recent NMR and neutron 
scattering experiments which reported anisotropic spin fluctuations. While CeCoIn 5 is a 
very complex system, we believe that the measurements of the field induced anisotropy 
provide a strong evidence for the symmetry of the superconducting gap. 

7. 2. k-(BED T- TTF) 2 Cu(NCS) 2 

The nature of the superconductivity in quasi-2D k-(BEDT-TTF) 2 X salts [in the 
following abbreviated as k-(ET) 2 X] , where the ion X can, for example, be Cu(SCN) 2 , 
Cu[N(CN) 2 ]Br or I 3 , has attracted considerable attention. In these compounds, the 
large molecules are coupled, forming narrow bands with low carrier density. It is known 
that ET molecules constitute a two dimensional conducting sheet in the crystal foe- 
plane, alternating with insulating layers of anions X. Within the conducting layer, ET 
molecules are arranged in dimer pairs with alternating orientations (Fig. 26). In k- 
(ET) 2 Cu(NCS) 2 , superconductivity occurs in proximity to the AF ordered state in the 
phase diagram, implying that the AF spin-fluctuations should play an important role 
for the occurrence of superconductivity; some (but not all) of the electronic properties 
of these superconductors are strikingly similar to the high-T c cuprates [146[ 1147] . 

The structure of the superconducting order parameter of k-(ET) 2 X salts has been 
examined by several techniques [Hj. Results strongly favoring ci-wave pairing with line 
nodes came from NMR [14811149] , thermal conductivity, and penetration depth [151(1152] 
experiments pSfl on X=Cu[N(CN) 2 ]Br and Cu(NCS) 2 . The STM [153] and mm- 
wave transmission [154] experiments reported strong modulation of the gap structure, 
although they arrived at very different conclusions regarding the nodal directions. In 
contrast to these experiments, specific heat measurements on K-(ET) 2 Cu[N(CN) 2 ]Br 
near T c suggested a full gap[155j. 

Figure [27] depicts the T-dependence of k. Since the phonon thermal conductivity, 
K ph , dominates the electronic contribution, n el , near T c , the enhancement of k below 
T c reflects the increase of the phonon mean free path as the electron pairs condense. 
Figures [28] (a) and (b) depict the if- dependence of k in perpendicular (HI. be— plane) 
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and parallel (H\\ be— plane) field, respectively. In the perpendicular field, k(H) shows a 
monotonic decrease up to H c2 above 1.6 K, which can be attributed to the suppression of 
the phonon mean free path by the introduction of the vortices [1501 1156] . Below 1.6 K, 
k(H) exhibits a dip below H c2 . The minimum of k(H) arises from the competition 
between n ph , which always decreases with H, and K el , which increases steeply near H c2 . 
Consequently the magnitude of the increase of k(H) below H c2 provides a lower limit of 
the electronic contribution, which grows rapidly below 0.7 K; k^/ K n is roughly estimated 
to be >5% at 0.7 K and >15% at 0.42 K, where and Kyi 3X6 the electronic and total 
thermal conductivity in the normal state above H c2 , respectively. 

We now move on to the angular variation of k as H is rotated within the 2D 
ki-plane. Figures [29] (a)-(c) display k(H, (f>) as a function of (ft =(q, H) at low 
temperatures. Above 0.72 K k(H, <fi) shows a minimum at = 90°, indicating simply 
that the transport is better for the heat current parallel to the vortices. On the other 
hand, at lower temperatures, the angular variation changes dramatically, exhibiting a 
double minimum as shown in FigsEH] (b) and (c). In all data, we fit as a sum of 
three terms: a constant, a two fold, k 2 0, and a fourfold, n^. These fits are shown by 
the solid lines in FigsEH] (a)-(c). Since a large twofold symmetry is observed even above 
0.7 K, where K ph dominates, k 2( ^ is mainly phononic in origin. In what follows, we will 
address the fourfold symmetry which is directly related to the electronic properties. 

Figures [291 (d)-(f) display K44 normalized by K n [66]. At T=0.72 K, the fourfold 
component is small: \C^\/K n < 0.1%. On the other hand, a clear fourfold component 
with \C^\/K n ~ 0.2% is resolved at 0.52 and 0.43 K. As discussed before, the 
contribution of n el grows rapidly below 0.7 K and contitutes a substantial portion of the 
total k at 0.4 K. Therefore it is natural to consider that the fourfold oscillation is purely 




T (K) 



Figure 27. Temperature dependence of the thermal conductivity in zero field. The 
heat current q was applied along the &-axis. Upper inset: The resistive transition at 
T c . Lower inset: The Fermi surface of k-(ET)2Cu(NCS)2- The Fermi surface consists 
of quasi-lD and 2D hole pocket. The node directions determined in our experiment 
are also shown. 
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Figure 28. Field dependence of the in-plane thermal conductivity (a) in perpendicular 
and (b) in parallel field [H || c) at low temperatures. Deviation from the horizontal 
line shown by arrows marks H C 2 . Inset: Field dependence of the electronic thermal 
conductivity An el in parallel field. 

electronic in origin. Although \C^\ at 0.42 K is as small as 0.2% in K n , it is likely close 
to 1.5-2% of and is also a few per cent of n el (0), assuming k^/k u ~ 0.15. The band 
structure of the crystal is very unlikely to be an origin of the fourfold symmetry, as the 
Fermi surface of this material is nearly elliptic with twofold symmetry, and the fourfold 
modulation is negligible, if present at all [157] . Hence we conclude that the observed 
fourfold symmetry originates from the superconducting gap nodes. 

The main question is whether the anisotropy of n el is associated with the DOS 
oscillation or Andreev scattering off the vortices. From general arguments [50, [53], the 
density of state effects dominate at low T, such as our experimental temperature of 
T c /30. Also, since K el increases with H, as shown in the inset of Fig. 2(b), we believe 
that the DOS enhancement underlies the if-dependence of K el at 0.42 K. 

If the DOS oscillations indeed dominate, K el attains its maximum value when H is 
directed along the antinodal directions, and has a minimum when H is along the nodes. 
According to Refs. [53| llllj . \C^\ in the d-wave superconductors arising from the DOS 
oscillation is roughly estimated to be a few percent of n el (0), which is in the same 
order to the experimental results. Since exhibits a maximum when H is applied 
parallel to the b- and c-axes of the crystal, we conclude that the gap nodes are along 
the directions rotated 45° relative to the b- and c-axes; the nodes are situated near the 
band gap between the ID and 2D bands (see the upper inset of Fig. 1). This result is 
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Figure 29. (a)-(c)Angular variation of K,(H,<j>) in \hqH\=2 T for different 
temperatures. cj> is the angle between q and H. The solid lines represent the result of 
the fitting by the function k(H, (f>) — Cq + C24, cos 20 + cos 40, where Co, C24, and 
C40 are constants, (d)-(f) The fourfold symmetry obtained from (a)-(c). 



consistent with the STM experiments [153J. 

We emphasize that the determined nodal structure is inconsistent with the recent 
theories based on the AF spin fluctuations. In the AF spin fluctuation scenario, it 
is natural to expect the nodes to be along the b- and c- directions since the AF 
ordering vector is along the b-axis, which provides partial nesting. If we take the 
same conventions for the magnetic Brillouin zone as in the high-T c cuprates with 
d x 2_ y 2 symmetry (see Fig. 1 (c) in Ref. [147j ). the superconducting gap symmetry of 
k-(ET) 2 Cu(NCS) 2 is d xy . Recently, it has been suggested that the nodal structure 
depends on the hopping integral between ET molecules, even if the superconductivity 
is mediated by AF fluctuation. For instance, d xy symmetry dominates over d x 2_ y 2 
when the dimerization of ET molecules is not too strong, which appears to be the 
case for k-(ET) 2 Cu(NCS) 2 |158j . Moreover, when the second nearest neighbor hopping 
integral tf, between the dimer is comparable to the nearest neighbor hopping integral 
t c (h ~ t c ), the (i^-symmetry is stabilized [159]. Indeed, in k-(ET) 2 Cu(NCS) 2 , tb/t c 
is 0.8, which is close to unity. It has also been argued that the charge fluctuations 
rather than spin fluctuations may be relevant to the unconventional superconductivity 
in k-(ET) 2 Cu(NCS) 2 . Consequently, our results for the gap symmetry should serve as 
a constraint on future development of the theories of superconductivity in this family 
of compounds. 
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7.3. Sr 2 Ru0 4 

Ever since its discovery in 1994 [8J, the superconducting properties of the layered 
ruthenate Sr 2 Ru04 attracted considerable interest |160| . The crystal structure of this 
material is same as La 2 Cu0 4 , a parent compound of high-T c cuprates (Fig. 30). 

The superconducting state of Sr 2 Ru04 stimulated great interest because NMR 
Knight shift remains unchanged from the normal state value below T c , indicating that 
the pairing state may be a spin triplet [21]. Recent phase sensitive experiments are 
controversial: although odd-parity of the superconducting wave function was suggested 
by Ref. [161j . it was also claimed that the result of Ref. [161] can be interpreted in 
the even parity framework [162j . /xSR measurements report the appearance of static 
spontaneous magnetic field below T c , which can be interpreted as a sign of broken time 
reversal symmetry [32]. The specific heat C r [163j . NMR relaxation rate |164] and thermal 
conductivity [165] indicate the presence of nodal lines in the superconducting gap. These 
results have motivated theorists to propose new models for the superconductivity in the 
ruthenates [166111671 1168U169] . We, of course, address this issue here from the standpoint 
of the thermal conductivity measurements. 

Inset of Figl3l7a) shows the T-dependence of k/T in zero field for crystals with 
T c =1.35 K and 1.5 K. The system is very pure, as is clear from the electrical resistivity of 
the order of 0.1fiQ-cm. At low T the electronic contribution to the thermal conductivity 
is dominant. At the superconducting transition, n/T(H) shows a kink. Figures |3~T1 (a) 
and (b) show the if -dependence of k for the sample with T c =1.45 K in perpendicular 
(H_L afe-plane) and parallel fields (H|| a6-plane), respectively. In both orientations, k 
increases with H after an initial decrease at low fields. The subsequent minimum is 
much less pronounced at lower temperatures. At low T, k increases linearly with H. 
For the in-plane field k rises very rapidly as H approaches H c2 and attains its normal 




Figure 30. Crystal structure of S^RuOzi. 
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Figure 31. Field dependence of the in-plane thermal conductivity (q|| [110] ) of 
Sr2Ru04 (T c =1.45 K) in (a)perpendicular field HI. afr-plane and (b)parallel field 
ff||[110]. In perpendicular field, k is indepedent of H below the lower critical fields. 
Inset: Temperature dependence of k/T in zero field for two crystals with different T c 
(T c =1.45 and 1.32 K). 



value with a large slope (dn/dH), while k in perpendicular field remains linear in H 
up to H c2 . In superconductors the slope of k(H) below H c2 increases with purity of 
the sample [156] , so that the data for the in-plane field suggest a clean limit. A rough 
estimate can be done as follows: let T be the pair breaking parameter estimated from 
the Abrikosov-Gorkov equation ^(1/2 + T/2irT c ) - ^(1/2) = ln(T c0 /T c ), where * is 
a digamma function and T c0 is the transition temperature in the absence of the pair 
breaking. Assuming T c0 =1.50 K and A = 1.76T C , T/A is estimated to be 0.025 (0.067) 
for T c =1.45 K(T C =1.37 K). Thus the dependence of k(H) observed in very clean crystals 
is consistent with the k of superconductors with line nodes. 

We now discuss the angular variation of the thermal conductivity. Figures [32] (a) 
and (b) depict k(H, <fi) as a function of = (q, H) [65]. In all the data k(H,4>) can 
decomposed as in Eq.(15). Figures [33] (a)-(d) show K^/Kn as a function of after the 
subtraction of Kq- and /^-term from k. Figure [341 depicts the if-dependence of \C^\. 
In the vicinity of H C 2 where k increases steeply, \C^\/k u is of the order of a several 
percent (see Fig. I3"37a)). We point out that both the sign and amplitude of in the 
vicinity of H c2 is mainly due to the in-plane anisotropy of H c2 . In FiglHH we plot the 
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Figure 32. (a) In-plane thermal conductivity normalized by the normal state value in 
H rotated within the afo-plane as a function of <j> of S^RuCU (T c =1.45 K). q is applied 
to the [110]-direction. (b) The same plot for the sample with T c =1.37 K. q is applied 
to the [100]-direction. <f> is the angle between H and q. The solid lines represent the 
twofold component in fc(0)//t n . 



amplitude of the fourfold oscillation calculated from the in-plane anisotropy of H c2 . The 
calculation reproduces the data. 
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Figure 33. (a)-(d) The fourfold symmetry n^/Kn at several fields obtained from 
Fig. [35] (a). 
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Figure 34. The amplitude of four fold symmetry as a function of H/H C 2- The 
filled circles and squares indicate \Cig\/K n at T=0.42 K and 0.55 K, respectively. The 
solid line represents \C\e\jn n calculated from the 4- fold symmetry of H C 2. Inset: In- 
dependence of k. The arrows indicates the points we measured C'^g. 



(CVI/^n decreases rapidly and is about 0.2-0.3% at lower field where k increases 
linearly with H (see Figs. [33] (b) and (c)). At very low field where k decreases with H, 
no discernible 4- fold oscillation is observed within the resolution of \C^\/ K n < 0.1% (see 
Fig. l33T d)). Thus the amplitude of the fourfold oscillation at low field of Sr 2 Ru0 4 is less 
than 1/20 of those in CeCoIn 5 and k-(ET) 2 Cu(NCS) 2 . These results lead us conclude 
that the line nodes are not located perpendicular to the plane, but located parallel ot the 
plane, i.e. horizontal node. If we accept the spin triplet superconducticity with broken 
time reversal symmetry, the gap symmetry which is most consistent with the in-plane 
variation of thermal conductivity is d(k) = Aoz(k x + ik y ) (cos ck z + a), in which the 
substantial portion of the Cooper pairs occurs between the neighboring Ru0 2 planes. 
Similar conclusion was obtained from the c-axis thermal conductivity measurements 
[170] . Recently, the in plane- variation of heat capacity has been reported down to 
100 niK. Below 200 mK, small but finite four-fold oscillation, which disappears at low 
H , was observed [72], [73] . This indicates that the gap function has a modulation around 
the c-axis, though finite gap remains. 

At an early stage, the gap symmetry of Sr 2 Ru0 4 was discussed in analogy with 3 He, 
where the Cooper pairs are formed by the exchange of ferromagnetic spin fluctuations 
[160J. However the inelastic neutron scattering experiments showed the existence of 
strong incommensurate antiferromagnetic correlations, and no sizable ferromagnetic spin 
fluctuations, indicating that the origin of the triplet pairing is not simply a ferromagnetic 
interaction [171] . The present results impose strong constraints on models that attempt 
to explain the mechanism of the triplet superconductivity. We finally comment on 
the orbital-dependent superconductivity scenario, in which three different bands have 
different superconducting gaps [172j . In this case, our main conclusion can be applicable 
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CeColn 5 K -(BEDT-TTF) 2 Cu(NCS) 2 Sr 2 Ru0 4 




YNi 2 B 2 C PrOs 4 Sb 12 UPd 2 AI 3 

Figure 35. The nodal structure of several unconventional superconductors, including 
quasi-two dimensional heavy fermion CeCoIns, organic k-(BEDT-TTF)2Cu(NCS)2, 
Ruthenate Si^RuO^The gap structure takes into account additional fourfold 
modulation as indicated by the specific heat measurements.), borocarbide YNI2B2C, 
heavy fermion PrOs4Sbi2, and heavy fermion UPd2Al3, which are determined by 
angular variation of the thermal conductivity 

to the band with the largest gap (presumably the 7-band). 
8. Summary 

In this paper, we discussed "how to determine the superconducting gap structure in 
the bulk?" . We show that the measurements thermal conductivity and specific heat in 
magnetic field rotating various directions relative to the crystal axis can determine 
the position and type of nodes. In Fig. [35], we summarize the nodal structure of 
several unconventional superconductors, including borocarbide YNi 2 B 2 C [67J, heavy 
fermions UPd 2 Al3|62j, CeCoIn 5 [6T], and PrOs4Sbi2[68j, organic superconductor, k- 
(BEDT-TTF) 2 Cu(NCS) 2 [66], and ruthenate Sr 2 Ru0 4 [65], which are determined by 
the present technique. While more theoretical work is clearly needed to arrive at a 
more quantitative description of the data, we feel confident that the method provides 
a uniquely powerful route towards determination of the nodal structure in novel 
superconductors. 
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